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SVtiOPSIS 

of the Ph. D. -thesis on 

A MICROCONTINUUM APPROACH TO SOME BIO-FLUID DYNAMICAL PROBLEMS 

D. PHILIP 

Department of Mathematics 
I.I.T. KANPUR 

In recent years study of human systems through interdisci- 
plinary collaboration has led to the fascinating area of 
Bio-fluid dynamics. This subject seeks to understand the 
mechanism of living systems using the concepts of fluid dynamics. 
Since physiological system in itself is complex, several 
simplifications are made to analyse the problem. One such 
simplification is to treat bio-fluid as Newtonian fluid. In most 
of the work till to-date bio-fluid dynamical problems have been 
analyzed using classical Navier-Stokes equation. But 
physiological fluids in general behave like suspensions, e.g. 
blood is a suspension of red cells, white cells and platelets in 
plasma; cervical mucus is a suspension of macromol ecul es in a 
water like liquid etc. Various theories, e.g. microcontinuum 
fluid theory [ Er ingen (1964), Ariman et al . (1974)], continuum 

theory of mixtures [Atkin and Craine (1976)], etc., have been 
suggested to account for the suspension behaviour of fluids. This 
has motivated the work in this thesis where an attempt has been 
made to analyse some physiological 
microcontinuum approach. ^ 


flow problems 


using 
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The following physiological flow situations have been 
analyzed in the thesis : 

(i) flow through mild stenosis, 

(ii) dispersion of a solute in a bio-fluid, 

(iii) peristaltic transport, and 

(iv) propulsion of spermatozoa in Cervical Canal. 

Ue have studied the above mentioned problems by considering 
relevant bio-fluid to be of simple microfluid [Kang and Eringen 
(1976)] nature which accounts for micro deformation and micro 
rotation of the micro particles (suspended particles). In the 
last problem, we have also taken micropolar fluid model for 
mucus . 

This thesis consists of 5 chapters. 

Chapter I is of introductory nature and provides motivation 
for the present work done in this thesis. A survey of the 
pertinent literature is given so that the importance of the 
problems analysed in the subsequent chapters can be seen in right 
perspective. It also includes the basic equations governing the 
flow of a simple microfluid. 

In Chapter II, we study the problem of flow through a 
stenosed artery. Here incompr esaible simple microfluid model for 
blood is assumed which account for the micro-def ormat ion and 
micro-rotation of the particles. The flow is assumed to be steady 
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and laminar. The analysis has been made for very mild stenosis 
and the effect of cell free layer near the tube wall has been 
accounted in the model. The effects of various parameters on flow 
resistance (impedence) and wall shear stress have been discussed. 
It is noted that the increase in the resistance due to the 
stenosis height gets enhanced for certain combinations of simple 
microfluid parameters. 

In Chapter III, the problem of dispersion of a chemically 
reactive solute in a simple microfluid is analysed considering 
the walls to be catalytic. The combined effects of homogeneous 
and heterogeneous chemical reaction is studied under Taylors 
limiting condition for flow through (i) channel and (ii) 
axi-symmetric tube. The effects of various parameters on 
equivalent dispersion coefficient have been discussed for both 
the cases. 

Chapter IV deals with the peristaltic motion of simple 
fflicrofluid through a distensible duct. The study consists of two 
parts (i) flow through a circular tube (ii) flow through a 
channel. In part I the analysis has been made for very long 
wavelength and hence the intertial terms in the equation of 
motion, governing the simple microfluid flow, h^ve been 
neglected. In part II, low Reynolds number flow of a simple 
microfluid through a distensible channel has been considered, 
under the long wavelength approximation. The regular perturbation 
analysis is used and the expression for stream function upto 
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first order has been obtained. The variation of pressure drop and 
other flow charact er ist i ca with respect to the simple microfluid 
parameters have been discussed. 

In Chapter V, we discuss the sel f -propulsion of spermatozoa 
through mucus filling the cervical canal , which is taken as a 
channel with flexible boundaries. The chapter consists of two 
parts : in part I, we have taken micropolar fluid model for mucus 
while in part II simple microfluid model has been considered. The 
spermatozoa is modelled as a two dimensional sheet which while 
swimming, sends down lateral waves of finite amplitude along its 
length. The model also considers the motion of the walls due to 
mascular activity which is represented by the transverse waves 
travelling along the flexible walls of the channel in the 
direction opposite the motion of the sheet. The analysis has been 
carried out for the inertia free flow under the assumption that 
the waves travelling along the channel walls and along the sheet 
are in synchronization under steady state and thus have same 
wavelength. The analysis presented here suggests that if the 
transverse waves are induced on the channel walls by some 
bio-chemical means, then it may be possible to reduce the speed 
of spermatozoa. This observation might be useful in fertility 
control . 

In our study here, the effect of simple microfluid has been 
discussed by taking arbitrary and independent variations of the 
microcont inuum viscosity coefficients, since no experimental data 
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could be ascribed to the simple microfluid parameters. Hence, 
the results presented here must be used with caution. It has been 
observed that though the qualitative behaviour of the flow 
characteristics for a given simple microfluid is similar to those 
for Newtonian fluid, they show quantitative variations. The 
results show, reasonably well, the capabilities of microcontinuum 
approach in the physiological flow problems. 
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CHAPTER I 


GENERAL INTRODUCTION 

"If you believe that a vital field of science is a source of 
important mathematical problem, then the field of bio-science is 
the proper domain for you. The field of biology is so rich and 
complex that one can not visualize its being exhausted in the 
next one or two hundred years.” 

Richard Bellman 


1.1 INTRCa>UCTION 

In the last few decades attempts to understand the human 
system through interdisciplinary collaboration have led to the 
fascinating area of Bio-fluid dynamics. Bio-fluid dynamics seeks 
to understand the mechanism of living systems by using the 
concepts of fluid mechanics. Since the physiological systems, in 
general, are very complex, theoretical analyses have been made 
using several simplifying assumptions. In most of the work which 
has appeared to-date the relevant fluid la assumed to be 
Newtonian and problems have been treated using classical 
Navier-Stokes equations. However, it is observed that most of the 
bio-fluids behave as suspensions of deformable or rigid particles 
in a Newtonian fluid, e.g. blood is a macroscopically continuous 
suspension of particulate matters such as red cells, white cells 
and platelets in a microscopically continuous fluid (blood 
plasma) [(Houssay, 1955)]; Synovial fluid is a suspension of 


hyaluronic acid molecules in plasma [Chandra (1980)]; 


the 
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cervicAl mucus is a suspension of raacromol ecul es in a water like 
liquid [Odeblad (1959, 1962)] etc. 

In view of this some researchers have tried to account for 
the suspension behaviour of the bio-fluids by considering 
non-Newtonian models while still using the classical continuum 
approach. However, for a suspension, the particles contained in 
the fluid, may rotate under the velocity gradient of the flow 
field; thus, in addition to the linear momentum, a spin angular 
momentum is also Introduced. This gives rise to the stress tensor 
becoming asymmetric, and hence the breakdown of the continuum 
hypothesis. Horeover, the classical approach becomes inadequate 
when the response of the body is sought to an external physical 
effect in which the length scale is comparable to the average 
molecular size. This limitation of the classical continuum 
approach has led to the development of theories of microcontinuum 
in which continuum media are now regarded as a set of structured 
particles which possess not only mass and velocity but also a 
substructure with which is associated a moment of inertia density 
and a micro-deformation tensor [Eringen (1964), Allen et al . 
(1967)]. These theories have been used to describe the mechanics 
of complex fluids like polymeric suspension, animal blood etc. 
and has found wide application in physiological and engineering 
problems [Cowin (1972), Arlman et al . (1974a), Kang and Eringen 
(1976)]. Therefore, in this thesis we have attempted to study the 
following physiological problems using the microcontinuum 


approach 
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1. Flow through a stenosed artery, 

2. Dispersion of a solute in a fluid, 

3. Peristaltic transport in a duct, 

4. Propulsion of spernatozoa in cervical canal. 

In our study here, we have taken simple microfluid model 
[Eringen (1964)] for the bio-fluid in each of the problem. In 
addition to it the last problem has also been analysed by 
considering the motion of the micro-organism in micropolar fluid 
flow, through a channel. 

Ue now give a brief introduction for simple microfluid and a 
brief survey of the literature pertaining to the above topics, so 
that the work done in this thesis can be seen in the proper 
perspective . 

l.a SEMPLE MICRCa^LUID 

The fundamental concept of classical continuum mechanics is 
that all material bodies possess continuous mass densities, and 
the laws of motion and the axioms of constitution are valid for 
each part of the bodies regardless of its size. Thus, such an 
approach considers media as a dense collection of point masses 
which are void of internal structures. The theory becomes 
inadequate when more refined and complete description of material 
is sought. This necessitates the considerations of the intrinsic 
motions of material constituent and has led to the development of 
the theory of microcontinuum in recent years [Eringen and Suhdbi 
(1964), Allen and Kline (1969), Kaloni and DeSilva (1969)]. A 
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review of various theories in this direction is given by Arinan 
et al . (1973a). The theory of microfluids, as introduced by 
Eringen (1964), deals with a class of fluids which exhibit 
certain microscopic effects arising from the local structure and 
micro motions of the fluid elements. These fluids can support 
stress moments and body moments and are influenced by the spin 
inertia. Thus, in the microstructured fluid model, material 
points of the body are considered to be small deformable bodies 
approximated by points with a triad of deformable directors 
attached to it. 

In the theory of simple microfluids while many of the 
principles of classical continuum mechanics remain valid, they 
have to be augumented with additional balance laws and 
constitutive relations. For example, the presence of microscopic 
elements in a fluid gives rise not only to classical Cauchy 
stresses but also to couple stresses due to micro element 
interactions; also the orientation of the micro element as well 
as the internal angular momentum due to their rotational motion 
needs to be considered. These balance laws and constitutive 
relations have been derived by Eringen (1964). We reproduce, 
below, the balance equations and the constitutive relations 
governing the flow of a simple microfluid [Eringen (1964)]. 
(These equations are given here in rectangular coordinate 
system). 

Conservation of mass : 

5 ? * = 


0 


( 1 . 1 ) 
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Balance of monentum 


‘kl.k * P “ 0 


( 1 . 2 ) 


Balance of first stress moiBents : 




Cl. 33 


Conversion of micro-inertia 


SI 


km 


-j*? ^ V - i u - i u - 0 

km , r r rm rk kr rai 


(1.4) 


The constitutive equations for incompressible simple microfluid 
are 


"^kl = -^^kl^^Aidj^l ^l^(kl) ^^tkl3 " ^kl^ 


(1.5) 


®kl ~ ^‘^kl ■*■ ^^‘^kl * ^^2^Ckl) 


(1.6) 


X, , = (y. V + y„ V ) <5, - + (y . u. + y^ v , ) <5, 

klm ^^1 mr,r ^2 rm.r'^ kl 4 lr,r ^5 rl,r^ km 


+ (^7 '‘^kr.r * ^8 ’^rkjr^ *^lm * ^10 ^kl,m * ^11 ^km,l 


* »'l2 "lk.« * ''13 %k.l * ''H “iB.k * ''15 %l.k 


where 


p = mass density, 

Vj^ = velocity vector. 


i, ■ = microinertia tensor, 
km 

t, , = stress tensor, 
kl 


fj^ = body force per unit mass, “ spin inertia. 
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stress av6ra£e tensor, 

X = the first stress laosient tensor, 
klm 


gyration tensor. 


^ Im ~ first body isoaent per unit mass, 


^ SS — f Y" 

kl 2 k,l 


''l.k’ ■ “kl “ 2 (''k.l * ’l.k’ 


( 1 - 8 ) 


and to , fJ, C,. Co *nd k are viscosity coefficients. Here 

partial di f f er ent iat ion is denoted by (,) and material time 
derivative by a superposed (.)• 


The symmetric part of the gyration tensor, represents 

iiiicro-def ormat ion while the skew symmetric part, ^[Icl] 
micro-rotation of the microelements . 


In our study in the subsequent chapters, we restrict our 
attention to two-dimensional motion of incompressible simple 
microfluid where the body forces and body couples are neglected. 

__ f 

Thus, the velocity V, fiyration tensor pressure P, have the 

forms 


= (U (X.Z) . 0 , 


W (X.Z)) 


-kiCX.Z) = 


^^(13) "^[13] 

when k=l. 

1 = 3 

» f 

'^(13) " "^[13] 

when k=3. 

1 = 1 

0 

otherwise . 



PCX,Z) 
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Thus the balance equations for incofflpressibl e fluid flow 
reduce to 


a? 




= P 


u + 


t a V 

C31) ^ 


SZ 




[31] 


az 


1 &^u 1 

2 axaz 2 TTI 

&L, 




u 

SZ 


Cl. 9) 


SF 

SZ 


+ P 


SX^ SZ^ 


= P 


K.S V> 


u + 




(13) ^ 




[13] 


SX 


1 S^U 1 s^u 

2 dXdZ 2 ^^2 

&K 


SU , 

^ " 51 " 


U ^ 


( 1 . 10 ) 


^ CV31) 




7 * 

S (v^3) 




...i * (K,-K3)— ♦ (Kj« 3)-— ^ . (Kj.Kp— ^ 


SX 


SX 


SZ 


SZ 


- 2(C,-2C3)v^^ 3^ t 2k 


' - 1 £H 1 £H 

^[13] 2 dZ 2 4»X 


= 0 


( 1 . 11 ) 


? ’ 

^ (- 31 ^ 


7 * 2 * 

s (V ) (V ) 

. y rjr TJ' 4 .-# 


^ ^- 31 ^ 




^X 


SX' 


SZ 


SZ 


* 2(Ci-2C2)-c 31) ^ 2k 


[31] 


1 ^ + i £M = 0 

2 ^X 2 dz 


( 1 . 12 ) 


^ 4 £H = 0 

SX SZ 


(1.13) 
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where 


■'l ■ * >'13 * >"15 * >'i * ’'12 * >'14 ^ ® 


*=2 * >'2 * >"11 * >"14 * >'5 ^ >'10 * >"15 ^ ° 


■'3 ' >'l * >'l3 * >'15 - (’’t * >"12 ^ >'14> ^ ° 


>=4 ' >'2 * >'ll ^ >'l4 ■ '>'5 * >"10 * >'l5> “ 


\ (1.14) 


To make an order of magnitude analysis, we nondimensional ize 
the various quantities as follows : 

(u.w) = (U,U)/V* , w = lly v! , p = Ph/(pV*\) 

iJ V 

X = X/L , z = Z/h, t = V*t/L, R* = phV*/p, 

K 

(E^,E2,E3) = C-Cj, ^2-^)/^, Kj = ^ h^ * h/L 


t * 

(where V is characteristic velocity, L is characteristic 

length along the axis of the channel, h is height of the 
channel ) . 


Thus, we get 


49x 


+ U 


*2 S 


&x 


a-. 




(13) 


az 


+ E. 




[13] 


az 


* *2 
£ & \J 

2 


1 a 


SxSz 


Sz 


* * 
R £ 


du ^ * du 

u -s— + £ W 
d'x Sz 


(1.15) 


Me 
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dx 3 




(13) _ 1 

Sx ~ 2 ax&z 2 ^ ^ 2 

Sx 


t t 

= R ^ 



w 


Sw 

Wz 


( 1 . 16 ) 


2 2 2 2 
- - *2 ^ ^31 - - *2 ^ ^13 - - ^ ^13 - - ^^31 

<^h-V‘ -rr- * -r¥ * " ('' 2 *'‘ 4 )— 4 ^ 


dx 


Sx 


»z 


dz 


ZCE^+ZE^) ^(13) + E 3 


^^[13] Jz ^ 


= 0 


(1.17) 


- - *2 ^^^^31 - - *2 ^^^^13 - - '*^^13 - - '^^’^31 

CK14K3). __ ^ CK^^K^)^ . CK^-K3)-^ 4 ck^-K,)-^ 


^x 


^x‘ 


Sz 


az‘ 


2(Ei42E2) + E 3 


. * Sn 

'^^[13] ~ ^ Jx * Jz 


= 0 


(1.18) 


du ^ * Sw 

Sx ^ Wz 


0 


(1.19) 


* 

Here it can be noted that for the inertia free flow (R =0) 

right hand side terms in eqns. (1.14) and (1.15) can be dropped. 

Further, for the case vhen h << L (i.e. infinite characteristic 

* 

length) terms containing £ can be dropped. 

1.2.1 Mlcropolar Fluid 


A subclass of simple microfluid which exhibits 
micro-inertial effects and micro-rotational Inertia is 


only the 
known as 


fflicropolar fluid. In micropolar fluids the directors attached to 
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the points are taken to be rigid which implies that the 
fflicrostructure have only micro-rotations and no microstretchings 
and shearings. For the sake of completeness we give below the 
constitutive equations for an incompressible micropolar fluid as 
given by Eringen (1966). 


^ij ‘ ■‘^ij ^ '''i.j " " * ''j.i - * ‘ijk “k 


i J If J Jf 1 


( 1 . 21 ) 


where and are stress tensor and couple stress tensor 

respectively, 6., is the Kronecker delta, an alternating 

1 j i jk 

tensor, p the isotropic pressure and the micro-rotation 

vector, p is the viscosity coefficient of the classical 
mechanics and x Is the new viscosity coefficient for the 
micropolar fluid. Comma after a suffix denotes partial 
di ££ er ent iat ion . 

The balance equations governing the flow of incompr essibl e 
micropolar fluid are [Eringen, 1966]: 

Conservation of mass : 


ST * ''k’-k ' ° 


Cl. 22) 


Conservation of linear momentum 


(P+*/2) - p,^ * p(C^ - V ‘ 0 
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Coneervation of angular moatentum ; 


* ‘’klm * ^'^k - ^ "^k' 


where is the body force per unit mass, 1 
couple and j is the nicro inertia constant, 
indicates the material derivative. 


= 0 (1.24) 

is the body 
Superposed dot 


1.2.2 Boxmdary Ck>ndltlons 


The determination of proper boundary conditions for 
microstretch and microrotation is a complicated topic, and a well 
posed problem requires a proof of uniqueness theorem which is not 
available even in the case of fluid dynamics [Kang and Eringen 
(1976)]. While the no slip conditions are assumed for the 
velocity components at the boundary surfaces, various types of 
conditions have been suggested for the microrotation [Migun 
(1984), Kirwan (1986)]. Thus, for the skew symmetric part of the 
gyration tensor one can assume : 


(a) perfect adherence at the rigid wall i.e. v., ,, = 0 

[kl] 

(b) moment free condition [Eringen (1966)] 

i.e. X-,,- n =0at the wall 
[kl]m ra 

(c) given fixed gyration at the wall [Ariman (1971)] 

i.e. = V , a constant 

[kl] o' 

(d) ~ ®''^kl 0 ^ a ^ 1 [Aero et al . (1965)] 


However, not much has been said about the micro-deformation 


of particles at the boundary surfaces. It seems that the 
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Bel ect ion 

of 

suitable 

boundary conditions 

for the 

simpl e 

iiicrofltiid 

depends upon 

the 

type of 

fluid 

and the 

physical 

aituat ion. 

In our work. 

we have aasuiaed 

the 

following 

type of 

condit ions 

as 

suggested by 

[ Kang & 

Eringen (1976)] 

for a 

auapsnsion . 









(13) 

= Vl3 

and 

"i> = 

[133 

*3“l3 


(1.25) 


A particular case o£ = A^ = 0 refers to a perfect adherence 
condition and denote that the fluid solid interaction is so 
strong that it does not allow the rotation and deformation of the 
microstructure at the wall. It represents the extreme case in 
which the walls act to completely inhibit microrotation and 
mi cro-def ormat ion . 

1.3 STENOSIS 

The word 'STENOSIS' means narrowing of the lumen of a tube 
or orifice. Thus, an arterial stenosis refers to a narrowed or 
constricted segment of an artery. An arterial stenosis is also 
referred to as 'atherosclerosis'. This can cause serious 
circulatory disorders by reducing or occluding the blood supply 
to an organ; e-g- , stenosis in coronary artery can cause 
myocardial infarction leading to heart failure, likewise in 
carotid artery it can bring about cerebral strokes. 

Though the process of initiation of arterial stenosis is not 
completely known, it is well recognised that many factors are 
involved in the development of stenosis, e.g. hypertension. 
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intake of sore lipid rich diet, lack of oxygen and lack of 
physical exercise [Spain (1966)]. It has been suggested that 
hydrodynamic factors such as wall shear stress, pressure 
distribution, separation phenomena etc. play an important role in 
the initiation and early progression of the stenosis [Rodkiewicz 
(1974), Stehbens (1982), Narem & Levesque (1987)] . Thus, in the 
recent years researchers have shown enormous interest to study 
blood flow through stenotic region. One of the earliest 
experimental study reported in the literature is attributed to 
Hann et al. [1938], who studied blood flow in a decreasing lumen 
of an artery. Young [1968] was the first one to study the problem 
mathematically. He considered a simplified linear model of flow 
in a mildly constricted tube, with the following geometry of the 
stenosis 


r‘cz, = . 

o 


^ {l - cos ^ (Z 
o o 

otherwise 



(1.26) 


where R(Z) is the radius of the artery in the stenotic region, R^ 

is the constant radius elsewhere, is the length of the 

stenosis and 6 is its maximum height. Stenosis location is 
a 

indicated by ,d, and Z the axial coordinate. Since then much 
attention has been given to the investigations of the 
steady /unsteady flow through arterial stenosis both analytically 
as well as experimentally [Forrester & Young (1970), Lee and Fung 
(1970), Caro et al . (1971), Cheng et al . (1972, 1973), Morgan and 
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Young (1974), Roth et al. (1976), Daly (1976), Deshpande (1976), 
Young (1979), McDonald (1979), Shukla et al . (1980a, b), Uilliams 
and Javadpour (1980) and Chaturani and Swamy (1983)]. Shukla et 
al . (1980c) have considered power-law and casson models for blood 
and have concluded that this non-Newtonian behaviour of the blood 
is helpful in the functioning of diseased arterial circulation. 
In recent years various micro-continuum models have been 
suggested to consider the suspension behaviour of blood. Hence, 
Sinha & Singh (1984), Srivastava (1985) have analyzed the effect 
of couple stress on the blood flow through a thin artery. They 
have shown that an increase in the couple stress parameter 
increases the resistance to the flow and the wall shear stress. 
Parvathamma and Devanathan (1983) have considered the steady flow 
of micropolar fluid through stenosed tube. They have shown that 
for a given Reynolds number, the shearing stress reaches a 
maximum immediately preceding the throat and then rapidly 
decreases in the diverging section. They have also shown that as 
the Reynolds number increases, the point of maximum shear stress 
moves upstream due to the slope of the tube wall. Recently 
Pralhad and Schultz (1988) have studied the stenosis problem for 
different diseases, considering blood as a polar fluid. Recently 
Sagayamary and Devanathan (1989) have studied the flow of couple 
stress fluid through tubes of varying cross-sections and have 
applied the results to the blood flow. 

Further it has been noted [Hershey et al . (1966), Middleman 
(1972), Lih (1975)], that when blood flows through arteries, near 
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the wall of the arteries, a cell free layer is formed i.e. there 
exist a plasma layer near the wall and the cells are radially 
distributed. Shukla et al . (1980a, b) have studied the effect of 
peripheral layer on the flow of blood to take into account the 
clustering of cells in the central zone and have concluded that 
such a layer helps the flow in a diseased system. 

Keeping in mind the above mentioned points, we consider in 
Chapter II, the blood flow through stenotic artery using micro- 
continuum approach. Blood is considered to be simple microfluid 
as suggested by Kang and Eringen [1976]. The effect of cell free 
layer near the wall is also considered. 

1.4 DISPERSI<»I IN BICft.X)GICAL FLUID FLOWS 

The process of dispersion is important in chemical as well 
as in biological systems. In particular, in all living beings 
matters like nutrients and metabolic products are transporated by 
blood as a result of diffusive and convective mechanisms; drug 
la transported as a result of diffusion through blood. Also the 
process of dispersion of a solute matter in a flow system is an 
important phenomenon in rivers, in oil pipelines, in hydruallc 
industrial devices and various chemical systems. Thus, various 
investigators have studied dispersion/diffusion under steady and 
unsteady conditions in different physical systems. 

Taylor (1953) was the first one to mathematically formulate 
the problem of dispersion. He used his model to measure the rate 
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of flow of blood in the arteries of animals by injecting a small 

volume of highly conducting liquid at some point. In his 

classical paper, Taylor (1953) discussed the dispersion of a 

soluble matter in a viscous fluid flowing in a ciricular pipe 

under laminar conditions. He found that relative to a plane 

moving with the mean speed of the flow, the solute is dispersed 

2-2 - 

with an apparent diffusion coefficient E v^/48D (where E, v^ and 
D are the radius of the pipe, the average velocity and the 
molecular diffusion coefficient respectively). In a subsequent 
study, Taylor (1954) showed that dispersion of solute can be used 
to measure molecular diffusion provided : 

4L 

F~ 

(L being a given length along the flow direction). Aris (1956) 
while extending the Taylor's work, found that the rate of growth 
of variance of the solute distribution la proportional to the sum 
of the molecular diffusion and Taylor diffusion coefficient. 

In the early studies it was assumed that the solute is 
chemically non-reactive. However, in a variety of problems in 
chemical engineering, diffusion of solute takes place in the 
presence of an irreversible first order chemical reaction, e.g. 
hydrolysis of eater, gas absorption in an agitated tank with 
chemical reaction [Bird et al . (I960)]. It is also suggested that 
the walls of the tube may be catalytic and hence may give rise to 
heterogeneous chemical reaction at the surface e.g., tubular flow 
reactors with heterogeneous catalysts. In view of this, Cleland 
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and Uilhelm (1956) have investigated the problem of a finite 
first order homogeneous reaction in a pipe under laminar flow 
conditions by a finite difference method and supproted their 
theoretical results with experimental data. Katz [1959] 
considered the chemical reaction catalyzed on the wall of the 
tube and discussed the heterogeneous reaction taking place at the 
wall on the concentration profiles. Later many generalizations/ 
extensions of these studies have been attempted to study 
dispersion of a solute in Newtonian and non-Newtonian fluid flows 
[Gumming et al . (1966), Gill et al . (1967, 1972), Ghosal (1971), 
Scherer et al . (1972), Gupta & Gupta (1972), Dutta et al . (1974), 
Shukla et al . (1979), Shukla & Parihar (19S0d), Chandra & Agarwal 
(1983), Das & Mazumder (1989) etc.]. 

Gill et al . in a series of papers (1967, 1970, 1972) have 
developed an exact method of analysis of unsteady convective 
diffusion system. Taylor's approach has been used to study the 
dispersion of a solute in non-Newtonian fluids by Fan and Hwang 
(1965), Fan and Wang (1966), Ghosal (1971), Shah and Cox (1974) 
for non reactive solute, while Dutta et al . (1974), and Shukla et 
al . (1979) considered dispersion of a solute in the presence of a 
homogeneous chemical reaction. Shukla et al . (1980d, 1982a) have 
studied the effect of peripheral layer on the Taylor dispersion 
coefficient and have found that the Taylor coefficient decreases 
as the viscosity or peripheral layer fluid decreases. Combined 
first order heterogeneous and homogeneous reactions were studied 
by Ufalker (1961), Solomon 4 Hudson (1967), Gupta 4 Gupta (1972) 
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and others. Gupta and Gupta (1972) in particular have considered 
unsteady dispersion of a solute with simultaneous chemical 
reaction in a liquid flowing in a channel under laminar 
conditions using Taylor's approach. Barton (1984) has used the 
Aris method of momenta to study the dispersion of reactive 
contaminant in pipe Poisuille flow in the presence of first order 
reactions at the wall. 

Using microcontinuum approach, Soundalgekar (1970) has 
analysed dispersion of a solute in the laminar fully developed 
flow of micropolar fluids flowing in a channel. In a 
subsequent paper Soundalgekar and Chaturnai (1980), studied 
the effect of couple stresses on the dispersion of a soluble 
matter in a pipe flow. Chandra and Agarwal (1983) studied the 
dispersion of a solute matter in simple microfluids flowing 
through channel and pipe under Taylor's limiting condition. They 
have considered the first order chemical reaction in fluid but 
did not consider the walls to be catalytic. Hence, in Chapter 
III, the dispersion of a solute in a simple microfluid is 
analysed considering the walls to be catalytic and the combined 
effect of homogeneous and heterogeneous reactions is considered. 

1.5 PERISTTALSIS 

The word peristalsis is derived from the Greek word 
"Peristaltikos" , which means clasping and compressing. It may be 
described as the process of fluid being transported through a 
distensible tube as contraction and expansion waves propagate 
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along its length on the boundary walls. Its laain function is to 
propel fluid contents in the tube. It plays an indispensibl e role 
in transporting many physiological fluids inside living bodies, 
e.g. urine transport froa kidney to bladder through ureter, in 
bile duct, ffioveioent of ovum in fallopian tubes, transport of 
spermatozoa in ductus efferentes of the male reproductive tract 
and in cervical canal of the female reproductive system, 
vasofflotion of small blood vessels, mixing the contents inside the 
gastro-intestinal tract etc. Many bio-mechanical and engineering 
devices have also been designed on the basis of the principle of 
peristaltic pumping. 

From fluid mechanical point of view, peristaltic pumping is 
characterized by the dynamic interaction of fluid flow with 
movement of flexible boundaries. Hence several authors have 
studied peristaltic transport of fluids from fluid dynamical as 
well as physiological points of view. Excellent reviews of the 
literature in this direction have been given by Jaffrin and 
Shapiro (1971), Rath (1980), Srivastava & Srivastava (1984). 

In particular Burns and Parkes (1967) have discussed the 
peristaltic motion in both two dimensional and axisymmetric cases 
under the assumption of small Reynolds number and linearized 
boundary conditions with perturbation technique. Barton and 
Raynor (1968) studied peristaltic motion in a circular tube with 
reference to small intestine using long as well as short 
wavelength approximations. Fung and Yih (1968) considered an 
idealized two dimensional analogy of the real problem and studied 
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the problem of peristaltic waves induced by a sinusoidal 
travelling wave motion of the walls at moderate amplitudes. Later 
this was extended by Yin and Fung (1969) to axisymmetric flow in 
a circular tube by perturbation method of solutions in terms of 
the ratio of the wave amplitude to tube radius. They had also 
taken into consideration the nonlinear convective terms. Chow 
(1970) obtained the solution of the peristaltic transport in a 
circular cylindrical tube in the form of power aeries expansion 
when the Reynolds number is small. Mittra and Prasad (1974) 
studied the interaction of Poiseuilie flow with peristaltic 
motion. Gupta and Seshadri (1976) have analysed the peristaltic 
transport in a non-uniform tube. Lykoudis and Roos (1970) studied 
the fluid mechanics of the ureter, from a lubrication theory 
point of view. Liron (1976) has presented a complete solution for 
peristaltic flow in a pipe and in a channel by means of an 
expansion in the Reynolds number and wave number for arbitrary 
wave shape. Gopalan and Devanathan (1983) have studied the 
pressure variation in a peristaltic motion with arbitrary wave 
shape when a catheter is inserted. 

In all these studies the effect of viscosity variation of 
the fluid has not been taken into account, which may play an 
important role in some physiological situations. To highlight 
this effect of peripheral layer, Shukla et al . (1980e) studied 
peristaltic transport of a Newtonian fluid with viscosity varying 
across the tube, under long wavelength approximation. The 
analysis has been applied and compared with the observed flow 



21 


rates in intestine and ductus efferentes of male reproductive 
tract, and the importance of the peripheral layer viscosity has 
been pointed out. 

Though it has been found experimentally that biofluids (e.g. 
human feces, bile, mucus etc) behave rheologically as 
non-Newtonian fluid [Odebald (1959, 1962), Han and Barnett 
(1973), Patel et al . (1973), Rodkiewicz et al. (1979)], only a 
few researchers have paid attention to study the peristaltic 
transport of non-Newtonian fluid. Raju and Devanathan (1972) 
studied the peristaltic flow of a power-low fluid. They have 
obtained power series solution for the stream function and the 
velocity components in terms of the amplitudes of deformation. 
They have also discussed the influence of the applied pressure 
gradient along with the non-Newtonian parameter on the stream 
lines. Picologlou et al . (1973) studied the peristaltic motion of 
human faeces characterized by a power-law model from the point 
of view of pressure generation in the contracting part and 
concluded that the peristalsis, in colon, though rare, is a major 
propulsive mechanism. Mittra and Prasad (1974) extended the work 
of Fung and Yih (1968) to study the effect of elastic or 
visco-elast ic wall on peristaltic flow. They have found the 
existence of "mean flow reversal”, both at the centre and at the 
boundaries. Raju and Devanathan (1974) extended their previous 
work to include the discussions about the peristaltic motion of 
the visco-elastic fluid. Radhakrishnamacharya (1982) has analysed 
peristaltic transport of a power-law fluid. In view of the 
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power-law behaviour of the chyme and the presence of peripheral 
mucus layer in intestines, Shukla et al . C1982b) studied the 
peristaltic transport of a power-law fluid whose consistency 
varies along the radial direction and investigated the effects of 
the peripheral layer consistency and the pseudoplastic nature of 
the bio-fluid on the flow characteristics. Kaimal (1978), 
Srivastava & Srivastava (1989) have studied the peristaltic 
motion of a suspension of rigid particles in a tube of arbitrary 
wave shape, at low Reynolds number under long wavelength 
approximation. An asymptotic method was developed for the 
solution of peristaltic transport of a viscous fluid with solid 
particles in a flexible tube of arbitrary cross-section by Shen 
et al . (1981). Biihme and Friedrich (1983) considered the 
mechanism of peristaltic transport of an incompressible 
visco-elastic fluid by means of an infinite train of sinusoidal 
waves travelling along the wall of the duet, in the case of a 
plane flow. Peristaltic flow of power-law non-Newtonian fluid 
containing spherical particles which are neutrally buoyant, is 
studied by Rath (1983). Devi and Devanathan (1975) have analyzed 
peristaltic motion of micropolar fluid under long wavelength 
approximation. Srivastava (1986) studied the peristaltic 
transport of a couple stress fluid under zero Reynolds number and 
long wavelength approximation. He compared the results with 
those for Newtonian fluid model to show that the magnitude of 
pressure rise under a given set of conditions is greater in the 
case of couple stress fluid. Therefore, in Chapter IV we consider 
the flow of simple microfluid in a circular tube and in a channel 
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so as to account for the mi croatretchings , shearings and 
micro-rotations . 

1.6 SWIMMII.R3 OF MI CRO>C»:6ANI SMS 

The study of self-propelling micro-organisms in liquids has 

received Considerable attention in the recent years due to its 

physiological applications such as motion of spermatozoa in 

cervical mucus. Spermatozoa are the tiny living cells (having 
-3 -5 

length 5x10 cm, diameter 10 cm) which have their own 
motility. They move towards the oviduct through the mucus, 
filling the cervical canal, by sending waves of lateral 
displacement down the tails. Its propulsion resembles to that of 
fish and its direction of movement is opposite to that of the 
propagation of the waves. 

Since the pioneering work of Taylor (1951) who studied the 
propulsion of spermatozoa by modelling it as a two dimensional 
infinite extensible sheet with sinusoidal wave travelling down 
its length, several researchers have analyzed the motion of self 
propelling micro-organisms [Hancock (1953), Reynolds (1965), 
Lighthill (1974), Shack and Lardner (1974), Shen et al . (1975), 

Pirronneau (1975), Rubinow (1975), Shukla et al . (1978), Sinha et 

al. (1982), Parihar (1983), Blake et al . (1983) Shukla et al . 

(1988)]. 

In particular Gray and Hancock (1955) obtained the velocity 
of propulsion of sea urchin spermatozoa with remarkable accuracy 
using simple dynamical arguments. The earlier studies considered 
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movement of isolated sheets in unbounded domain. Lately some 
studies have been made to analyze the hydrodynamics of a swimming 
spermatozoa in a bounded fluid which may actually be the case 
once the sperm enters the cervical canal. Reynolds (1965) 
considered two dimensional model and used Taylor’s sheet approach 
to consider the effect of a rigid wall adjacent to the sheet. 
This analysis is based on the assumption that the amplitude of 
the tail motion is small. He observed that the transport of sperm 
through the cervical mucus is too rapid to be accounted for by 
the invitro swimming speed of approximately 50 pm/sec of sperm in 
cervical mucus. Llghthill (1975) presented a general theory of 
swimming of micro-organisms at low Reynolds numbers. Shukla et 
al . (1978) studied the swimming of spermatozoa in cervix by 
considering the interaction of the micelles aligned along the 
wall. They have also pointed out that the propelling velocity of 
the sperm increases due to the dynamical interaction of the 
cervix and due to the presence of the peripheral layer. Attempts 
by other reaearehera were alao made to explain the motion of 
spermatozoa in the female genital tract by considering the 
dynamical interaction of the wall. Smelser et al . (1974) 
investigated the swimming of sperm in an active channel in which 
the walls of the channel also vibrate. Blake et al . (1983) 
presented a theoretical model of ovum transport in the oviduct 
incorporating transport mechanisms due to ciliary and mascular 
activity of the wall by adding a force distribution term in the 
equation of motion. The contraction and relaxation of smooth 
muscles of uterus and the oviduct may give rise to peristaltic 
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motion of the channel wall. This effect was considered by Shukla 
et al. (1988). 

There is a strong likelihood that the motion of spermatozoa 
through the cervical canal may be studied from a micro-continuum 
view point. It has been observed that cervical mucus is a 
suspension of high molecular weight macromolecules in a water 
like liquid with viscosity 0.03P; Odebald (1959, 1962). In luteal 
phase the mucus resembles like a close mesh, having a spacing of 
0.03 fjm in them, Elstein (1971), Davajan (1971). This motivated 
Sinha et al . (1982) to study the propulsion of spermatozoa in a 
micropolar fluid. They have considered the channel walla to be 
rigid. Thus in Chapter V, swimming of micro-organism is studied. 
Ue mathematically analyze the motion of the spermatozoa in the 
cervical canal by considering the transverse waves along its tail 
and peristaltic motions along the cervical wall. Cervical mucus 
is considered as micropolar fluid in Part I and as a simple 
microfluid in Part II. 

1 . 7 PRESEMT STUDY 

In view of the above, we present here a microcontinuum 
approach to study some physiological flow problems. This approach 
accounts for the microrotations and microdeformations of the 
suspended particles. 

In chapter II, we study the problem of flow through a 
stenosed artery. Here incompressible simple microfluid model for 
blood is assumed. The flow is assumed to be steady and laminar. 
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The analysis has been made for very mild stenosis and the effect 
of cell free layer near the tube wall has been accounted in the 
model. The effects of various parameters on flow resistance 
(impedance) and wall shear stress have been discussed. It is 
noted that the increase in the resistance due to the stenosis 
height gets enhanced for certain combinations of simple 
microfluid parameters. 

In chapter III, the problem of dispersion of a chemically 
reactive solute in a simple microfluid is analysed considering 
the walls to be catalytic. The combined effects of homogeneous 
and heterogeneous chemical reaction is studied under Taylor's 
limiting condition for flow through (i) channel and (ii) 
axi-symmetric tube. 

Chapter IV deals with the peristaltic motion of simple 
microfluid through a distensible duct. The study consists of two 
parts (i) flow through a circular tube (ii) flow through a 
channel. In part I, the analysis has been made for very long 
wavelength by neglecting Inertial terms. In part II, low Reynolds 
number flow of a simple microfluid through distensible channel 
has been considered, under the long wavelength approximation. The 
regular perturbation analysis is used and the expression for 
stream function upto first order has been obtained. 

In Chapter V, we discuss the self-propulsion of spermatozoa 
through mucus filling the cervical canal, which is taken as a 
channel with flexible boundaries. The chapter consists of two 
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parts : in part I , we have taken aicropolar fluid stodel for otucus 
while in part II aisaple microfluid model has been considered. The 
spermatozoa is modelled as a two dimensional sheet which while 
swimming, sends down lateral waves of finite amplitude along its 
length. The model also considers the motion of the walls due to 
mascular activity which is represented by the transverse waves 
travelling along the flexible walls of the channel in the 
direction opposite the motion of the sheet. The analysis has been 
carried out for the Inertia free flow under the assumption that 
the waves travelling along the channel walla and along the sheet 
are in synchronization under steady state and thus have same 
wavelength and wave speed. 

In our study here, the effect of simple microfluid has been 
discussed by taking arbitrary and independent variations of the 
microcontinuum viscosity coefficients, since no experimental data 


could be ascribed 

to the simple microfluid parameters. 

Hence 

, the 

results presented 

here 

must be used with 

caution. It 

has 

been 

observed that though 

the qualitative 

behaviour of 

the 

flow 


characteristics for a given simple microfluid is similar to those 
for Newtonian fluid, they show quantitative variations. The 
results show, reasonably well, the capabilities of microcontinuum 
approach in the physiological flow problems. 



CHAPTER II 


FLOW OF SIMPLE MICROFLUID THROUGH AN ARTERY WITH MILD STENOSIS 

2.1 INTRODUCTION 

It is well known now that flow through stenosed artery can 
cause serious circulatory disorders. Hence, the flow 

characteristics of blood through a stenosed artery, have been 
extensively investigated in recent years. One of the earlier 
attempts to theoretically analyse the interaction of stenosis 
with fluid dynamics of blood flow was made by Young (196S) in the 
case of mild stenosis. This has been subsequently investigated 
by many other workers, for steady and unsteady cases [Forrester 
and Young (1970), Lee and Fung (1970), Morgan and Young (1974), 
Shukla et al . (1979)]. However it may be noted that in these 
studies blood has been characterised as Newtonian fluid and 
little attention has been given to the suspension nature of 
blood. The experimental studies have indicated that, under 
certain flow conditions, rheology of blood significantly deviates 
from the Newtonian fluid behaviour due to its suspension nature 
[Charm & Kurland (1965), Uhitmore (1968)]. For example, non 
parabolic velocity profile, existense of peripheral layer etc., 
have been observed for the flow through tubes of small diameter 
[Bugliarello 4 Sevilla (1970), Cokelet (1972), Goldsmith (1975)]. 
Thus attempts have been made to Theologically describe blood as a 
suspension of neutrally buoyant deformable particles in viscous 
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fluid using microcontinuum approach [Kline et al . (1968), Cowin 
(1972), Popel et al . (1974), Ariman et al . (1971,1973b), Kang and 
Eringen ( 1976 ) ] . 

In view of this Sinha and Singh (1984), Srivastava (1985) 
considered couple stress fluid model to study blood flow through 
stenosed tube. Though their models account for the size effects 
of the particles, it does not account for the angular velocity 
and deformability of the suspended particles (Cells). Pralhad 
and Shultz (1988) have discussed blood flow through a stenosed 
artery for different diseases by considering polar fluid model 
for blood, which accounts for the particles’ spin. However, the 
effect of cell deformability, which has pronounced effect on the 
apparent viscosity, is not considered in this model. Recently, 
Kang and Eringen (1976) have proposed simple microfluid model for 
blood. This model accounts for the microspin and micro-stretch 
of the particles. Therefore in this chapter, an attempt is made 
to analyse blood flow through a narrow artery with mild stenosis 
by considering simple microfluid model for blood. The model 
presented here considers only the localized effects of the 
stenosis and accounts for cell free layer of blood plasma near 
the tube wall [Bugliarello and Sevilla (1970), Shukla et al . 
(1980a, b)]. Further we assume, flow to be steady, as the 
pulsatile behaviour of blood flow is not very significant in 


narrow arteries. 
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2. 2 FORMULATIC»l 


Ue consider here two layer — model for blood flow in atenosed 
artery. The model consists of a central layer of simple 
microfluid in the core region and a peripheral layer of Newtonian 
fluid (plasma) near the wall of the tube. The geometrical 
configuration of stenosis, which is not very well defined, is 
assumed to be axially symmetric and the radius ECz) (Fig. 2.1) of 
the stenosed tube is given as follows [Young (1968)] 

d < z < L +d 
o 


R - 
o 


1 2n . 

1 + cos 1 — (z 

o 




R(z) = i 


R otherwise 

o 


( 2 . 1 ) 


Here R(z) is the 

radius o£ 

the tube 

with stenosis. 

R 

o 

is 

the 

constant radius. 

L is the 
o 

length of 

the stenosis and 

S 

B 

is 

the 


maximum height of the stenosis. Ue further assume that : 

(i) Fluids in both the layers are incompressible, 

(ii) Body forces and body moments are absent. 


(iii) Flow, which is due to the pressure gradient dp/dz, is 
one-dimensional and axi-symmetr ic , 

(iv) Stenosis is very mild, which implies that the variation of 
all the fluid quantities (except pressure) along the axial 
direction is negligible, 

v) The length of the tube is large compared to its radius 
(which on order of magnitude analysis implies that the 
inertial terms can be neglected). 
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(vi) The peripheral layer is of constant thickness, S. Thus the 
core region is given by 0 < r ^ R^(z) = R(z) “ S and 
peripheral region is given by R^(z) r R(z). 

Using the above assumption we now write the euqations 
governing the steady laminar flow in the two regions : 

Ca3 CC^lE REGICW CO < r < R^Cz33 : As mentioned above, here we 

have a simple microfluid flowing under the pressure gradient 

dp/dz. The field quantitites are, velocity component, v = (0, 0, 

V ), gyration tensor components and (other components 

p 13 31 

being zero) and pressure, p. Thus the equation of balance of 
linear momentum as given by equation (1.2) in Chapter I give 


S 

Sr 


ct^^) 




( 2 . 2 ) 


while the equations due to balance of first fflooient stress (1-3) 
reduce to 
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where t. s. ., X. are £iven by the constitutive equations (1.5 

ijijijk 

to 1.7) in Chapter I. The non zero coaponents can be written as. 
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. “'13*‘'31 _ ., . , “13-*'31 (2.15) 
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Ci» Co» y-i to >'^c the viscosity co-efficients. C » 

1 z 1 10 1 

C- and k have the dimensions of viscosity fj , while y.’s are of 

Z X 

2 

the dimensions pL . These co-efficients satisfy the following 
restrictions [Kang and Eringen (1976)], 
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p > 0 , 2 C 2 -C 1 2: k > 0, 4p(C2 - j - CC^/2)^ 2: 0 

( 2 . 16 ) 

^ ^2 ^ ^15 - ^ ^11 ^ ^14 - ° • 


^4 defined in chapter I. 


Here and £ive the symmetric and skew symmetric parts 

CldJ [Idj 

of the ^j^2‘ symmetric part i.e. ^( 3 ^ 3 ) represents local 

stretchings and shearings, while the skew symmetric part, 
give local intrinsic angular velocity of the substructur e - 


Cb) 

PERIPHERAL 

REGION CR^CzI 

< r < RCz>:) ! 

The 

equation 

for 

the 

Newtonian 

fluid (plasma) 

flow with V 

Ii 

O 

P 

and 


constant viscosity /j under pressure gradient is. 


1 

S 


r 

Sr 

L J 



- 


(2.17) 


THE BOUNDARY CONDITIONS 

(i) No slip condition gives, 

v^ = 0 at >1 = R(z) . (2.18) 


(ii) Symmetry of the flow implies that. 


Sv 

c 

St 


"(13) " ° "[13) ' “ at r = 0 


(2.19) 
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(ill) Continuity of velocity at the interface implies. 


V (R, ) = V (R, ) 


C2.20) 


(iv) Continuity of shear stress at the interface gives 


^v 
I 

Sr 


r=R. 


= t 


13 


r=R. 


( 2 . 21 ) 


(v) The determination of proper boundary conditions for ^ c^-. and 

V Ki J 

V.,,, at the interface is rather difficult. 

[kl] 

Here, following Kang and Eringen (1976), we assume similar 
condition for the microatretching and microspin at the 
interface. Thus, at r = Rj^(2), we have 


C13) 


[13] 


= A 


= A 


^v 

c 

1 ^ 

tfv 

c 

2 Sr 


( 2 . 22 ) 


(2.23) 


where A^^ and A^ are constants. 


a. 3 ANALYSTS 


Integrating the equation (2.2) using the equation (2.5) and 
the boundary condition (2.19) at r = 0, we obtain 


13 


1 “(IS) * *' ‘■C13J 


rP 

2 


(2.24) 
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which gives the interface condition ( 2 . 21 ) for stress in the 
following form : 


Sv 
I 

dr 


r=R. 


= t 


13 


RlP 


r=R. 


(2.25) 


Equation (2.17) alongwith the boundary conditions (2.18) and 
(2.25) gives v^ as. 


^ (r2 

i(j 


r^) 


Rj r 3^ R 


(2.26) 


Now substituting for X 223 , t 3 ^. 83 ^. X^ 3 ^, X 232 , t ^3 and 

from eqns . (2.5) to (2.12) in eqns. (2.3) and (2.4), we get 


7^ ..3J . 7^ . 2(Ci-2C2) ^(31, 


1 

"2*= '“'[31) - 1-5?^ 


= 0 (2.27) 


(K3-K3) 72 .33 . (K3-K,) 7 ^ . 2((^ - 2C3) .^33^ 




where V s 


^ .Id _ 

3 . 2 r dr 2 
dr r 


On adding (2.27) and (2.28) we get 


•'l “31 *2 “13 * «1 - ^ 2 ^ tvj 3 * .- 33 ) 


0 


(2.29) 
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while subtraction of (2.28) from (2.27) gives 
2 

*"3 ’ -^31 * *4 ‘'13 - “'ll - ‘*31 * 


Using eqn.(2.24), the eqns . (2.29) and (2.30) can be rewritten as 
foil ows : 


>=l' - »1 3 -ai ‘ *^2 t - “2 3 ‘'13 = ■> 


(2.31) 


K3 [ V 


2 , r, r „2 2 , 

"3 J ’^31 ^ *^4 t ^ - ^4 3 ^13 


krP 


(2^1 + k) 


(2.32) 


where 


V 2 

Ki 


*^3 “3 


K 

4 4 


(C, - 2C2) 


„ 2 
*^2 “2 


k HiJ + C3^)/(2/j + k) ^ 


k HiJ - Ci)/(2p + k) 


(2.33) 


The equations (2.31) and (2.32) suggest the following forms for 
^13 ^31 


"13 = ^^1 ^1 ^ ^ 2^2 ^1 ^ ^ - 


(2.34) 


31 


^1 ^1 ^2 ^1 ~ 


(2.35) 


where 


r V ^-2 -2 . ,-?T .—2 -2 . 

^1 ~*^4 " «4)/K3(ot - 013) 


£2 = -^2 C/9^- «2)/K^((3^ - oc^) 


(2.36) 
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a. = a. R , a = a R , ft = ft R , K. = tz 

X i o o o 1 „2 

u R 

o 


(2.37) 


2 2 

and a and /3 ar® obtained usinc the equations (2.31) and (2.32) 
and are £iven as the positive roots of the quadratic equation 


(K^K^ - K^K^) + x[Ccx^ + oi^) - (oi^ + «^) K^K^] 


^ 22 „„ 22 „„ 

* “2 “3 ''2'‘3 - “l “‘4 ''l’‘4 


C2.3e) 


[Ij^Cx) is the modified Bessel fimction of order 1]. 

Subst itutinfi the expressions for and in eqn. (2.24) and 

integrating we get. 


''c ' *>3 * i ^ ^ 


+ / (1 . fj.) - I ^ (1 - I^COTJ 


^2 Pr^ 

" 2 - - ^ 2 > 


(2.39) 


where is constant of intecration. 


The constants b^, b^ and b^ are determined using conditions 
(2.19), (2.21) and (2.22) along with the equations (2.23) and 


(2.25). Thus, the field quantities ^(^ 3 ^ (micro-stretchings), 

V,..-, (micro-spins) and the velocity v are obtained in the 
I u j c 

following form : 
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PR 


(13) 


4p 

PR 


V 


[13] 


4^ 


bi(l + fi) Ij(otr) + b2Cl + f2) I 


biCl-£i) I^(oir) + b2(l-f2) I 


- rj(2. 


(2.40) 


41) 


Her 6 , 


PR , , PR R, 

4^ T 


I^(otR^) - I^C(otr) 


} 


+ d,^I^(^Rp - l^aftry 


(2.42) 


d. = 
1 


-b . e . 
1 1 

a.I^(aiRi) 


bi = |^(l + 2A2)g2 - 2Aj^h2j / - 63^1] 




b^ = [(l+2A2)g^ - 2A^h 


1 ] > [« 1^2 " « 2 ^l] 


g^ = (-1) 


= (- 1 ) 


i + 1 


(l+f^) 


i + 1 


(1-f .) 


.) 

■) 


£. 

1 


e . = 
1 




«^) + (a^ - 01^) K2K3 


K 1 K 3 («2 


-2. 

" 3 ^ 


(E 1 -E 3 ) + (£^+£ 3 ) f. 

( 2 +E 3 ) 


Y (2.43) 
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R.(z) 

^ ; R = RCz)/R^ 

o 

d = d/L ■ and L = L /L. 

o o 




I 


l + cos2n (z-d - -j") 


)■ 


It may be remarked that 6=0 i.e. 
of no peripheral layer and the 
microfluid, while R^(z) = 0, gives 
microfluid. 


= 

R(z) gives 

the case 

tube 

is full of 

simpl e 

no core 

region of 

simpl e 


The flow flux Q which is defined as. 


R 


0 = 2fz 


I 


r V 


dr + Zn 


1 


V dr 
P 


can be obtained in the following form using equations (2.25) and 
equation (2.42) 


fiR 

0 = ^ ^ 


° F [R(z), R^(2)] 


(2.44) 


where 


F(z) = F[R(z), R^(z)3 


R^ 


(2 + E2)|^*^i|^o^“^1^ 


2Ij(«R^)'| 

1 





i^C/?Ri) 


2I^(/?R^)^ 


(2.45) 
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Using equation (2.44) alongwith the conditions 


P = P, 


at z = 0 4 P = P at z = L, 

^ o 


we get, 


Po"Pi 


Is 

_ BjjQ r dz 

nR^ Jo F(z) 


(2.46) 


It may be noted that F = F(R , R - 6) is constant in the region 

o o 

o < z < d, d+L < z < L. 
o 


Thus , 


_ _ _ BfjQ 

P_ - Pj - T 

o 1 ^4 

fiR 

o 


d+L 


L-L r 

-- — 

F(R=R ) i 


dz 

F(z) 


d 


(2.47) 


The peripheral resistance X, is given in the non-dimensional 
form as [Young (1968)], 


X = 


1-L 


F(R=1) 


d+L 


1 


dz 

F(2) 


(2.48) 


where 


X = (p.-p^)/Q . X = 8pL/(nR ) 
1 o o o 


(2.49) 


The equation (2.48) under the transformation 

1 1 J Os 

- (2-d-^) 

L ^ 

o 
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can be written as 


1 



l-L 

+ r 

d^ 


F(R=1) - 

F(^) 


o 


(2.50) 


The shearing strehs at the wall r in nondisiensional for® is 

w 

given by 


T 

W 



4/jQ ^w 



(2.51) 


a. 4 RESULTS AND DlSCUSSIC^f 

The resistance to the flow, X and the wall shear stress, t 

w 

are the two important characteristics of blood flow through 

atenosed artery. Therefore in the following we discuss the 

effect of various parameters on these flow characterist ics . X 

and T , given by eqns . (2.50) and (2.51) r espect ively , are the 
w 

non-dimensional forma of X. and The integral appearing in 

eqn. (2.50) is not amenable to analytical form, so we numerically 

integrate it using Simpson’s rule. Polynomial approximations 

are taken for the Bessel functions I^(x), I^(x). Apart from the 

usual dimensionless parameters, i.e. S (the maximum height of 

a 

the stenosis), L (the length of the stenosis) and S (the 

o 

peripheraql layer thickness), we encounter simple microfluid 
parametrs (E^, , E^ , , K^) in the expresssions of 

Xand T . 

V 
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The parameters Ej^, and E^ correspond to the 

deformability and rigidity of the particles, e.g. at fixed (j 

higher values of k represent more rigid structure, while higher 

values of and |C 2 j represent more flexible structure. 

However the precise rates at which these parameters increase or 

decrease with deformability or rigidity are not known. Further, 

these constitutive co- ef f i ci ents depend upon the initial 

concentration of the particles as well as on the particle shape, 

deformability, and the fluid viscosity; but again their explicit 

relations are not known. In view of this, we study the effect of 

independent variation of these parameters on the flow 

characteristics. The parameters and are the linear 

combinations of the viscosity co-ef f ici ents to have 

2 

the dimensions of juL . Thus the non-dimensional parameters ( = 

2 - - 
K./^jL ) for a narrow artery have been chosen as = 1.1 K^= 1.0, 

= -0.9, = 0.6 [Kang and Eringen (1976)]. 

In table I and II we discuss the case when =0 i.e. 

the microspin and the microstretch at the interface are zero 

while figs. (2.2 to 2.11) are for A^^ = A^ = 0.25. It is clear 

from Table 1 and 2 that both X and t (calculated at the maximum 

w 

height of the stenosis) increase as 6^ increases, but decrease as 
the peripheral layer thickness increases. This is similar to the 
behaviour of these characteristics as observed by Shukla et al . 
(1979) who took Newtonian fluids in both the regions but with a 
different viscosity. The resistance to the flow gets enhanced 
with the increase in and E^ but with decrease in E^ - 
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Fig. (2.2 to 2.4) show the variation of X vs <5 when <5 

8 

(peripheral layer thickness) = 0.05. It is observed that in the 

case of siaple aicroflttid also, the resistance to the flow 

becomes more as S (height of the stenosis) or L (len£th of the 

& o 

stenosis) increases . 

Fig. (2.5) shows the effect of the increase in peripheral 
layer thickness for a given simple microfltiid. It is clear that 
increase in peripheral layer thickness decreases X. Further X is 
more when the tube is filled with simple microfluid only in 
comparison to the case of a Newtonian fluid filled tube. Thus 
the presence of peripheral layer helps blood flow in diseased 

cas es . 

The effect of the simple microfluid parameters and are 
opposite to each other. It is clear from the figures that an 
increase in E^ enhances the resistance to the flow while that in 
E^ reduces it. Unlike the case of = A 2 = 0 where the effect 

of E. and E« are not significant, here X shows appreciable 

X « 

variations with and E^- The effect of the parameter E^, in 
this case, is to reduce X, as E^ varies between 0.1 to 0.2 (refer 
Figs 2.3 and 2.4), however X increases for higher values of E^ 
(Figs. 2.3 ad 2.4). The combined effects of L and <5 with the 
simple microfluid parameters are further elaborated in figs. ( 2.6 
to 2.8). 

In figs. (2.9 to 2.11) the shear stress at the maximum 

height of the stenosis, is plotted for various values of E^ , E 2 
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FIG. 2. 2 VARIATION OF A WITH Ss/Rq FOR DIFFERENT 
El AND Eo. 






SO'O 
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FIG. 2. 7 VARIATION OF ^ WITH Sg/Ro FOR DIFFERENT E, 

Lq/l and S. 






0.05 
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FIG. 2. 10 VARIATION OF WITH 6-/R„ FOR DIFFERENT 
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E^. The qualitative effect of these parameters on is 

similar to what had been observed in the case of X. It is noted 

from fi£. (2.11) that t increases as 6 increases and this 

V s 

effect is further enhanced in the case of simple microfluid. 


Table I 


PERIPHERAL LAYER THICKNESS C<5) = 0.1 


6 

a 


^2 

^3 

X 

T 

W 

1 . 

0.04 

0.5 

0.3 

0.07 

1 .1386791 

1.2284889 

2. 

0.04 

0.5 

0.3 

0.09 

1 .1704005 

1 . 2609274 

3. 

0.04 

0.5 

0.5 

0.07 

1.386752 

1 . 2284844 

fl . 

0.04 

0.5 

0.5 

0.09 

1 .1703953 

1 . 2609215 

5. 

0.04 

1.0 

0.3 

0.07 

1.1388472 

1 .2286735 

6 . 

0.04 

1.0 

0.3 

0.09 

1.1706287 

1.2611773 

7. 

0.04 

1.0 

0.5 

0.07 

1.1388398 

1.2286652 

s. 

0.04 

1.0 

0.5 

0.09 

1 . 706185 

1.2611660 

9 . 

0.08 

0.5 

0.3 

0.07 

1.1917793 

1.3878380 

10. 

0.08 

0.5 

0.3 

0.09 

1.2242435 

1.4217867 

11. 

0.08 

0.5 

0.5 

0.07 

1.1917751 

1.3878330 

12. 

0.08 

0.5 

0.5 

0.09 

1.2242379 

1.4217803 

13. 

0.08 

1.0 

0.3 

0.07 

1.1919561 

1.3880483 

14 . 

0.08 

1.0 

0.3 

0.09 

1 . 2244831 

1.4220702 

15. 

0.08 

1.0 

0.5 

0.07 

1.1919482 

1.3880389 

16. 

0.08 

1.0 

0.5 

0.09 

1.2244724 

1 .4220572 
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Table II 





PERIPHERAL 

LAYER THICKNESS (^) = 0.2 



6 

a 


^2 

^3 

X 

T 

W 

1 . 

0.04 

0.5 

0.3 

0.07 

1.0936822 

1.1811299 

2 . 

0.04 

0.5 

0.3 

0.09 

1.1095525 

1 .1970510 

3. 

0.04 

0.5 

0.5 

0.07 

1.0936797 

1 . 1811270 

4. 

0.04 

0.5 

0.6 

0.09 

1.1096492 

1 . 1970474 

5. 

0.04 

1 . 0 

0.3 

0.07 

1.0937872 

1 . 1812405 

6 . 

0.04 

1 . 0 

0.3 

0.09 

1.1096913 

1 .1971972 

7 . 

0.04 

0.5 

0.3 

0.07 

1.0937823 

1 . 1812353 

a. 

0.04 

1.0 

0.5 

0.09 

1.1096849 

1 .1971902 

9 . 

0.08 

0.5 

0.3 

0-07 

1.1451915 

1.3362423 

10. 

O.OS 

0.5 

0.3 

0.09 

1.1613095 

1 . 3524337 

11 . 

0.08 

0.5 

0.5 

0.07 

1.1451889 

1.3362394 

12 . 

0.08 

0.5 

0.5 

0.09 

1.1613060 

1.3524296 

13. 

0.08 

1.0 

0.3 

0.07 

1.1453001 

1.3363623 

14 . 

O.OS 

1.0 

0.3 

0.09 

1 . 1614529 

1.3525913 

15. 

O.OS 

1.0 

0.5 

0.07 

1.1452950 

1.3363567 

16. 

0.08 

1.0 

0.5 

0.09 

1.1614462 

1.3525839 
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2. 5 CCM4CLUSION 

The effect of simple microfluid flow (in the core region) 

with & Newtonian fluid peripheral layer in a tube with very mild 

stenosis has been investigated. It is observed that resistance to 

the flow (X) and wall shear stress (t ) increase when <5 

w ® 

increases. These characterist ics are further enhanced by certain 
combinations of simple microfluid parameters. Hence it may be 
inferred that, this suspension nature does not help blood flow in 
the functioning of the diseased arterial system. 



CHAPTER - III 


EFFECTS OF HETEROGENEOUS MD HOMOGE^C:OUS REACTIONS ON THE 
DISPERSION OF A SOLUTE IN SIMPLE MICROFLUID 

3.1 INTK<»>UCTION 

The dispersion of a. solute in fluid flowing through 
channela/pipes is important in chemical as well as biological 
systems. In one of the early studies, Taylor (1953) presented a 
simple mathematical model to study dispersion of a solute through 
a fluid. He observed that, relative to a plane moving with the 
average speed of the flow, the solute disperses with an 
equivalent dispersion coefficient which depends upon (i) the 
average speed of the flow, (ii) the radius of the tube, and (iii) 
the molecular diffusion coefficient. In his analysis Taylor 
(1953) assumed that the solute does not chemically react with the 
fluid. However, subsequently many investigators analysed the 
dispersion problem by considering first order homogeneous 
reaction, under laminar flow conditions. Katz [1959] discussed 
the influence of the chemical reaction catalyzed on the wall of 
the tube. The combined effects of homogeneous and heterogeneous 
chemical reaction for a solute dispersing in Newtonian fluid 
flow have been discussed by Ualker £1961], Solomon and Hudson 
[1967], Gupta and Gupta [1972] and others. 



L 


(q) Channel 



(b) Tube 


F(G-3-1 GEOMETRY OF THE PROBLEM 
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Here the equation la written by neglecting the axial diffusion 
and the term ~Kc represents volume rate of disappearance of the 
solute- D, is constant molecular diffusion coefficient, K is the 
first order homogeneous chemical reaction rate constant and v is 
the velocity of the simple microfluid along the z-axis. 


Further, assuming that the walls of the channel are 
catalytic (hence, allow a heterogeneous chemical reaction) the 
differential material balance at the walls gives [Katz (1959), 
Gupta and Gupta (1972)], 


-D 


Sc 

3Fx 


fc 


at X = h 




(3.2) 


where, fc gives the surface reaction rate. 


For convenience the equation(3 . 1 ) is written relative to the 
axis moving with mean fluid velocity (v) (with co-ordinates 
z* = z - vt , X* = x) which on using the dimensionless parameters. 


becomes 


e 


tv 

TT 


z 


z-vt 

“L 



1 ^ v-v ^ 

t W L 


S^C y 

D - Kc 


(3.3) 


(3.4) 


where L is the characteristic length along the flow direction and 
V is average fluid velocity given by ^ 
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V 


h 


H j • 

o 


(3.5) 


AsaujBing that the Taylor'a limiting condition ia valid i.e. 
the partial equi 1 ibr iiiin is established in any croaa-sect ion of 
the channel , eqn. (3.4) reduces to 


S^c 2 

>T -me 


Syi 




(V - V) 


where — is independent of x, ^ = 0 and 

S~z dt D 


(3.6) 


(3.7) 


characterizes the chemical reaction rate. 


The boundary condition (3.2) can be rewritten as 


— + yc = 0 at x=l 


— — - j'c = 0 at X = -1 

^x 


1 


(3.8) 


where ^ 



the surface reaction rate constant. 


Now to solve the equation (3.6) for concentration ,c, we 
need the expression of velocity ,v, of the simple microfluid. As, 
the flow is laminar and one dimensional we take v = ( 0 , 0 , v(z)) 
and the non-zero components of the tensor are ^ 13 ( 2 ) 

v^^(z) only, thus the governing equations [Chapter 1, ( eqns . 1.9 

to 1 . 12)3 for the flow of simple microfluid in a channel can be 
written as [Chandra and Agarwal (1983)], 
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1 * v” 

T- v” - E V + E, + v' 

h 1 C13) 3 ^2h (13) 


) - hP = 0 


(3.9) 


^31 ^ - 2(Ei-2E2) -^, 3 ^ 


- ^ 3 (h- 2 -^^ 33 ) = 0 


(3.10) 


(K 1 -K 3 ) + (K^-K^) ^”3 2 (E^+ 2 E 2 ) ^(13) 


^ "3fir * 2 -[13]> = ° 


(3.11) 


Here ( ) = — K^, K 3 , are defined as in Chapter I. 

dx 


Here *^rin symmetric and skew-symmetric parts of 

f [.loj 

respect ively , as defined by equation (2.15). 

X o 

The boundary conditions for v, ^^ 3 ^^ are taken as follows 

[Kang and Elngen (1976), Chandra and Agarwal (1983)] 


(i) V = 0 ; = ir = *'[13] ■ ir *'’ ** ’ 


(ii) V = 0 ; ^(13) " ° 


*■[13] = “ 


at x =0 


(3.12) 


Solving the above equations (3.9 to 3.11) as in Chapter II, 
the expressions for ^£13} ^ obtained as 
[Chandra and Agarwal (1983)], 


(13) 


t 


Ph ''1 Sinhotx 

m 2~ 'l'’ 


Sinha 




inty^x I 

Sinh^ J 


(3.13) 
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(13) 




X + 


b 

2 1 Sinhc* 2 


^ Sinh^ 


(3.14) 


and 


V = - 


2aj 


- 2 * 

1-x +d^ (Coshot 


coshotx) + d 2 (cosh/? 


cosh/^x) 


(3.15) 


where 


* 

d. 

i 




[-2Aih2 + - e^h^J 


* * * * 


>2 = r-2Ajhj t (lt2A2)«^)/[gjh2 - g^h^J 


b . e . 

X 1 

a^sinha. 


= (-1) 


i + l 


i + 1 


(1+f*) 


- - ""l^i > 


= (-1) " { 


ci-f^ 

- ^2^i> 


(af - * (4 - ^2^ ^2^3 

^1^3 ^^1 - "3^ 


(E. - 


E3) . (E^ . E3) f. 


(2 + E 3 ) 


(3.16) 


and a^^ = a ; a.^ = ft and other constants such as ot^ , etc. are 


as given in Chapter II but non-dimensional ized by h instead of 


R . 
o 



65 


Thus the average velocity of the fluid is obtained as 


V = 


Ph‘ 

Zfj 


’’ 

" -1 


r — " 

7 

4 + d 

3 “l 

cosh^ - 

* 

^ ^2 

. - sinh/9 

cosh/3 - — 


^ . 


ft J. 


(3.17) 


Solving equation (3.6) with equations (3.15), (3.17) and the 

boundary conditions (3.12), we get 


fi 


\m m m 


1 -f 


6co3hrox 

^^(y .m) 


sinha . 


^ d* 

i^l ^ 


. a, coshinx 

1 i 


2 2 2 
m ni(a.-m )sinhinj 


cosha . X 


a ^ jn. ^ 

— fH 


a , - HI 

1 


★ 

c ^ 




2incoshinx (in,oi,/9) 


~ 2sinhro F^Cm,ot,ft) 


+ 2y sinhm F^(in,ot,/?) 


( 3 . 18 ) 


where 


a. hv . 

* o Sc 


D ’ "o = 


Ph‘ 


Fj^(in,oi,/3) 


F3(n. 


- 2 ) 2 : 
2 j £ — I a. -in 


2 


(3.19) 




* 

d. a. 


•-—I u. , o , 

I 

^ a . -in 


sinha. 
2 1 


(3.20) 


F3(in,oi,/3) = 


d .sinha . 


i_ ,1 - 1_ + V ^ ^ 

2 ^3 2 ) a 

HI m ^ — ' i 


(3.21) 



66 


= ^sinhiB {in^ + coshmj 


C3 . 22) 


[Note : Throughout this chapter the suinfflation index ,i, takes the 
value 1 and 2 ] . 

The average solute flux, Q, across the plane which moves 
with the mean speed of the flow can be obtained from, 


jk 

Q = 2h J c (x) (v-v)dx . 
o 


(3.23) 


Substituting the expression of v and v from eqn. (3.15) and eqn. 
(3.17) in eqn. (3.23) we get. 


n v^h ^ 

ih = >'> 


(3.24) 


where 


(m,oi,/3,?') = 


1 * Sinh a ^ * Sinh ft 

3 * 'll ^ * <*2 ^ 


a 


ft 


F - F (m,«,/^) + —3 

Dm o 


* - * — 
dj^Sinhot d^Sinh/S 


Sm'' oi(oi^-m^) /3(^^-m^) 


1 2 2 
F. a- Sinhm - ^ Coshm + Sinhm) 


5 "m 


m ni 


d*Sinha d* Sinh(5 

11 2 


F (m.oi,/ 3 ) + — - 27 ^,52 „2, , 

^ 3 m «(« -ffl ) "ffl ) J 
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2 . 


2 

i = l 


d. Sinha. 

F Cm, c*. ft) (-1— i) - 


pinh(a^+m) Sinh(a^-m)^ 


a , +01 
1 


(a; 


-o) J 


m 


— Sinha. - — Cosha. + — :r Sinha. 

a , 1^2 1 3 X 

1 a . a . 

i 1 


Tk * 


^ ‘^1*^2 fsinhCa+Z?) SinhCa-ft) ^ 1 ^ 

I « 


a 


OL+ft 


/- 2 2 , 

(a.- ft) a^-m J 


*2 

- Sinh 2., * I 

(a -m ) L 1 


C3.25) 


and Fg = p —y |^|2>'niFj^(0i,oi,/5) - 2ffiF2(Bi,a,/5)j- Coshm 

- ^ZyF^CmfCH,^) - 2r^F ^Cm,a,ft)'^ Sinhm 

Now cooiparing the equation (3.24) with the Fick’a 
diffusion, namely 


(3.26) 


law of 


J* = - D* ^ (3.27) 

oz 

* 

we get the equivalent dispersion coefficient D with which the 
solute disperses relative to a plane moving with the mean speed 
of the flow la obtained as follows 


* 

D 



M^(m,ol,/3,^) 


( 3 . 28 ) 
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3.3 PART II : FLOW THROUGH PIPES 


In this aection we consider the combined effects of 
homogeneous and heterogeneous chemical reaction on the dispersion 
of a solute in a simple microfluid flowing through a circular 
hollow cylindrical tube of radius R (fig- 3-lb). Ue shall analyse 
this problem under the same assumptions as discussed in Part I- 


Thus , the non- di mens iOjg^. form of concentration equation 


moving coordinates under Taylor's limiting condition can 
written as. 


in 

be 


1 


- Sc^ 2 -R 

Sc 



(r — ) -me = (» 

w 

r 

Sc 

Sc 

Sz 


(3.29) 


where , 


2 


m 


KR^ * 
D ’ ^ 


(z-vt)/R, 



independent of r and 



The boundary 

Sv 

Sc 

Sr 


onditions are 

0 at 

yc = 0 at 


[Katz (1959)], 



(3.30) 


where 

Rf 

For this one dimensional flow under the uniform pressure 
gradient 4^= P, the velocity, v of a simple microfluid is given 
by eqn. (2.42) [in Chapter II, with <5=0 and Rj^(z) = R(z) = R^]. 
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Thus 


PR 


2 ^ 


V = 


4aj 


(r^-1) 




Io(«)-I^Car) 




(3.31) 


The expressions for and other constants are as given in 

eqn . (2.43). The average velocity of the flow, v, is given by. 


J 


2nrv dr 


= V 


i • .. -¥;}•.{ 


^ ^2 - 


2I^(/5) 1 

ft 


(3.32) 


where 


v^ = PR“/4Ai 
o o 


From eqns . (3.31) and (3.32) we get. 


v-v = V 


^ - r +d 




(a) 


♦ dj. 




= v^ Gj^(r) (say). 


(3.33) 


Using eqn. (3.33) the equation governing the solute 
concentration can be written as. 



(3.34) 



Thus the aolution of the above eqn. (3.34) alongwith the boundary 
conditions (3.30) gives the concentration of the solute in the 
following form : 


c(r) = c 


u r, _ 


yl^(mr)' 

lm^[ 

2 - 2 

SI 2m 



(2+y) 

f (?^,ffi) 2 

6 HI 


JU 

- 1 


i = l 
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-1. 

■ 2Ij(ap 

' I^(nir)^ 


[■ailjCa.) * yl„(a.) 
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a . m 
*- 1 

F.(^',m) 
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2 2 
a, - m 
^ 1 


2.. 

I (a.r) 1 
+ 1 

2 a . 

m ^ J 

2 2 
a . - m 

1 aJ 


where 


Hr Hr 

c 


Rv 

o 

D 


Sc 

Hr 

Sz 


and 


F^(^ ,10) 


ml^(in) + j'I^(m) 


(3.35) 


The average solute flux Q across the plane which moves with the 
mean speed of the flow can be obtained from, 


1 


Q = 


J 


r (v - v) c(r) dr 


(3.36) 


by substituting the expression for v, v from equations (3.31) and 
(3.32) respectively. Thus we get. 
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wher e 


(roiy ,oi) + H.(ia,y,^3) 


* ‘‘l '^2 




+ d 


H ^ Cm.y , ot ) 


+ d : 


H ^ Cb . j ' ,/ 9 ) 


(3.37) 


where 
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4 y 2 + J ' 

“4 “T“ 

m in 


2 m 




SI 
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(3.38) 


1 f _ 41 («) I («) 
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ml j^(in)I^(a) - oilj^(oi) I^(b) 2 1 (a) I (m)' 
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( 8 y + 4 m ^ ^^ m ^) 

2 ffl^ F ^( r , m ) 


f2I^(a) Ij(ot) 4I^(a) 


.-2 2 . -2 
(« -m ) I ot 


2 oi 


-3 
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< Mlj ^( a ) + y IqC ") 

=:^2 2 
a - m 


1 r ^ 

IB^O J IB 


( ra ) Ij ( b ) 
2 2m 


4I^(ffl) 


F ^( r , m ) 


(3.39) 
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3.4 RESULTS AND DISCUSSION 

It may be pointed out here, that for the particular case 
when y —0 the equations (3.18) and (3.35) reduce to the case of 
Chandra and A^arwal (1983), while m=0 refers to the situation 
when the dispersion does not undergo any irreversible chemical 
reaction in the fluid, but has surface reaction at the wall. 

The effect of various parameters such as chemical reaction 

rate constant, m, surface reaction rate constant, y, and the 

viscosity coefficients on the equivalent dispersion 

coefficient can be seen through the function , in the case of 

channel flow and for the pipe flow. Hence these expressions 

have been numerically calculated using Simpson’s rule for 

different seta of values of various parameters. The results are 

presented by choosing the coefficients K^’s as = 1.1, = 

1.0, K. = -0.9 and K, = 0.6 and A = A. = 0.25. 

’3 4 12 

The effect of various parameters on are shown in Figs* 

(3.2 to 3.9) and on in Figs. (3.10 to 3.14). Fig. (3.2) shows 
that, for a given simple microfluid, decreases as the chemical 
reaction rate constant *m’ Increases. This effect is enhanced as 
the surface reaction rate constant increases. The effect of y on 
is further elaborated in Figs. (3.6 to 3.9). It is observed 
from these figures that shows slight decrease as r increases 
from 0 to 2 and then approaches asymptotic value as r is 
increased beyond 2. This increase in with y. becomes more 

appreciable as m Increases. In the case of the tube similar 




Aa usual comparing eqn. (3.37) with Pick's law of diffusion we 

* 

find the equivalent dispersion coefficient D with which the 

solute disperses relative to a plane moving with the mean 

* 

velocity of flow, D is given by 



(3.42) 
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FIG. 3.2 VARIATION OF WITH m FOR DIFFERENT 
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FIG.3.4 VARIATION OF Mi WITH m FOR DIFFERENT Ei 8. E 2 
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FIG.3.5 VARIATION OF M) WITH m FOR DIFFERENT 
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FIG. 3.9 VARIATION OF Mi WITH FOR DIFFERENT m , Ei & E 3 
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effects of surface reaction rate constant y on are observed 

(Figs . 3 . 10 to 3 .13) . 

The effects of siaple microf luid paraiaeters in the presence 
of surface reaction at the wall are depicted in Figs. (3.3 to 
3.5) (with y = 4.0). These figures show that increases as the 

viscosity coefficient increases, but decreases as the 

viscosity coefficient and E^ increase. However, variation of 
with E^ is not very significant and is almost negligible for 
smaller values of E^^. This behaviour is similar to the case of no 
surface reaction rate at the wall (y = 0) [Chandra and Agarwal 
(1983)]. Combined effects of the viscosity coefficients E^ , E^ , 
E^ and y for m = 0.5 can also be observed through Figs. (3.7 to 
3.9). In case of tube flow the variation of with respect to 

various parameter shows similar behaviour as has been observed in 
the case of channel flow. However, approaches asymptotic 

values for higher values of y (around y - 5). Also unlike the 

channel case here, the effect of on is not very 

insigni f leant . 

3. 5 CONCLUSICW 

The combined effect of homogeneous and heterogeneous 

reaction rate constants on D is to decrease it for a given set 

* 

of simple microf luid parameters. The decrease in D with 
homogeneous reaction rate becomes more in the presence of surface 
reaction on the walls. The equivalent dispersion coefficient 
decreases as the simple microfluid parameters E^ and E^ increases 
but it increases with the parameter E^ . 
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FIG. 3.11 VARIATION OF M 2 WITH FOR DIFFERENT 
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FIG. 3.12 VARIATION OF Mo WITH 75' FOR DIFFERENT 






chapter IV 


PERISTALTIC TRANSPORT OF A SIMPLE MICROFLUID 

4,. 1 IMTRODUCTKWI 

The study of peristaltic transport of a fluid has received 
considerable attention during last two decades because of its 
potential applications in biological systems. The earlier studies 
were primarily conducted with a view to understand the mechanism 
of urine transport from kidney to bladder. Therefore, several 
authors have analysed peristaltic motion of a Newtonian fluid, 
under various kinds of approximations [Burns and Parkes (1967), 
Barton and Raynor (1968), Fung and Yih (1969), Chow (1970), Li 
(1970), Jaffrin and Shapiro (1971), Mittra and Prasad (1974), 
Manton (1975)]. However, in many other physiological situations 
the physiological fluid behaves like suapenslon/non-Newtonian 
fluid. Hence in recent years attention has been given to study 
the peristaltic transport of a non-Newtonian fluid in 
channels /pipes . In particular, many studies have been carried out 
with reference to vasomotion of the small blood vessels, by 
considering blood as visco-elastic fluid [Raju and Devanathan 
(1972), Bohme and Friedrich (1983)], power law fluid [Raju and 
Devanathan (1972), Radhakrishnamacharya (1982), Shukla et al . 
(1983)]. Attempts were also made to analyse peristaltic transport 
of a suspension [Kaimal (1978), Shen et al. (1981), Rath (1983), 
Srivastava and Srivastava (1989)]. These studies do not consider 
the particle size effect or rotation and stretchings of the 
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particles. Renuka Ravlndran & Devi (1970), Devi and Devanathan 
(1975) analysed the problem for micropolar fluids to account for 
the microrotational effects of the particles, while Srivastava 
(1986) considered couple stress fluids to see the particle size 
effect on peristaltic motion. 

In this chapter we consider the peristaltic transport of a 
simple fflicrofluid with sinusoidal wave travelling down the walls 
of the duct. This chapter consists of two parts. In part I, we 
have analysed the peristaltic motion of a simple microfluid 
through a cylindrical tube by using very long wavelength 
approximation and neglecting the inertial terms [Shapiro et al . 
(1969), Shukla et al . (19S2b)]. In part II, we study the two 
dimensional flow of a simple microfluid through a channel under 
peristaltic wave. Approximate solution, upto the first order, is 
obtained by perturbation analysis with slope parameter as the 
perturbation parameter [Radhakrishnamacharya (1982)]. The effects 
of simple fflicrofluid parameters on the pressure rise across the 
duct have been discussed. 

4.2 Part I : FLOW THROUGH A CYLINDRICAL TUBE 

Ue consider here the peristaltic transport of a simple 
microfluid filled in a cylindrical tube. The walla of the tube 
are assumed to be executing sinusoidal wave motion due to 
peristalsis and so the geometry of the tube wall can be described 
as (Fig. 4.1), 
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R 



sin [ 


2n 


(z 


- Ct) 


( 4 . 1 ) 


where , is the radius of the tube in the absence of peristaltic 
wave, b is the amplitude of the wave, X the wavelen£th ,C, the 
wave propagation velocity, t is the time and z’ is the axial 
co-ordinat e . 


4. 2. 1 Fonrailation 

The cylindrical tube is filled with a homogeneous and 
incompressible simple microfluid and the flow is assumed to be 
axisymmet r 1 c . The constitutive equation and equation of motion 
are as given in Chapter I. In what follows, we consider a moving 
frame of reference whose velocity is synchronised with the 
peristatic wave speed ,C, and moving in the positive axial 

f 

direction. Thus transforming the stationary co-ordinates (r and 

» t > 

z ) to the moving co-ordinates (r,z), we get r = r, z = z+Ct . 

f t 

Further, if (u , o, v ) are the velocity coaponenta in the 
stationary co-ordinatea and (u, o, v) are those in moving 

f 9 

co-ordinates , then, u = u, v = v+C. 

Now, using very long wavelength approximation and 

neglecting inertial terms [Shapiro et al . (1969) , Shtikla et al . 
(1980e)3 equations of motion in the moving co-ordinates reduce to 
the following dimensionless form 


iE = 0 

9r 


(4.2) 
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(4.3) 


(4.4) 


K. 


* 

r It 1 

r 
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- “'I "31 * *4 
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7^ - 

“4 1 

' ^^13 
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• 




dp 


2 + E 


(4.5) 


3 dz 


where 


z -Ct 


t = 


Ct 


r = 


- 1 ; V 


(13) C *^(13)’ ‘"[13] 


R 

o 

c“ "^[133 


(4.6) 


p R . - 2 

P = rj^ ; R = ^ = 1 + £• Sin 2nz; '7 

^ o 


d" 1 d 1 

*-2 - ~ -2 
r #r r 


_ 

and -o = ^ < 1 is the amplitude ratio, K., E. as given in Chapter 

I and II respect ively. 


[For the sake of siiapllcity the bars over the flow variables have 
been omitted here.] 

The no-slip condition on velocity component gives, 

V = -1 at r = R(z) - (4.7) 

The boundary conditions on micro-stretch and micro-spin 
components are assumed as given by Kang and Eringen (1976), i.e.. 
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V 


(13) 


= A — and v 


=:A„— atr = R(z). (4.8) 


[13] 


2 


The axiaymmetry of the flow gives. 


dv 

Sr 


" ° ’ ’^(IS) " '"[IS] " ° 


at r = 0 


(4.9) 


is 


From equation (4.2), it is seen that the pressure p 
independent of r and it is only weekly dependent on t [Shukla et 

al . (1980)] i-e. 

- /-s (4.10) 

p = p( 2 , t) = p(z) 


In the following, we shall write. 


^ = P . 

Sz 


Th» ^qustlona (1.3 to l.B) aro aolv.d along with the 
bonndary oondltiona, which give the velocity v aa , 


.1 e I 




I (^R)-I^(/3r) 

O u 




(4.11) 


[Have d, and d, ace a. given in Chapter ii with 5, replaced hy i] 
The flux, q* in the .ovlng ayaten, ia given in di.enaionleaa for. 
by, 
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q 


nR^c 

O 


- 1 


2rv dr , 


(4.12) 


Uhile flux, Q , in the stationary frame is given by, 


R(z ,t) 


0 = 2n 


J 


f f f 

r V dr , 


(4.13) 


* 2 

which in dimensionless form Q (= Q /nR^C) can be written as 


R(2) 

Q = 2 1 r (v+1) dr 


J 


q+R' 


(4.14) 


Using (4.11) in (4.12) we get. 


-4 

-2 PR^ 
q = -R - -j- 




RI,(«R). 

^ ^ A 1 


1 \ 

1 * 

[ 2 - ^ 0 ^“^^ 

a J 

+ d- -I 

fR^ - 

(/5R) 

RI^(^RY 

2 1 

0 
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(4.15) 


which gives. 


£p = p = -i- 1 ) 


dz 


F(R) R" R 


(4.16) 
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where 


FCR) 





RIjCoiR) 


a. 


} 


+ d. 




RIjC^R) 

ft 



Further, the time averaged flux Q, (In dimenaionlesa form), for a 
complete time period T = X/c is obtained by. 


0 


J »di 

o 


(4.17) 


which on uaing equation (4.14) and the expression for R (equation 
4.6) gives , 

2 

Q = q+i + . (4.18) 

The pressure rise, 

Ap = p(0) - p(X) (4.19) 

across one wavelength is the same whether measured in the fixed 
or moving co-ordinate system, thus, it can be calculated from the 
equation (4.16) using the relation : 



1 



Ap 

= - j < al > ■ 

(4. 

20) 


o 



Further, it 

mqlY b6 pointed out that. 

in moving frame 

of 


reference, the pressure and flow appear stationary and therefore 
the flow rate ,q, measured in the moving co-ordinates is a 
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constant that varies neither with time nor with position along 
the tube axis [Shapiro et al . (1969)]. Hence, equation (4.16) and 
(4.20) give, 


where 


Ap = Ij ql^ 



(4.21) 


(4.22) 


and 



(4.23) 


Using equation (4.21) in equation (4.18) we get the expression 
for 0 in terms of Ap as, 


Q = 




(4.24) 


FRICTION FORCE s 

* 

The dimensionless friction force (= ; F being the 

friction force at the wall in the stationary co-ordinate system 
which is same as in the moving system) across one wavelength can 
be obtained as , 

1 

y = r (- ^ ) di 

J dz 

o 
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which on subaitutlon of equation (4.16) giv 


es 


5^=1:, + ql. 


where 


= J 


3 I F 
o 


dz 


(4.25) 


(4.26) 


4r. 2. 2 RESULTS AND DISCUSSICM4 

The effect of the flow parameters on pressure rise and 
friction force is seen through equations (4.21) and (4.25). These 
equations involve integrals I and which are not amenable 
to analytical solutions and hence have been numerically 
integrated using Simpson’s rule and polynomial approximations for 
the Bessel functions. The results are presented by taking Kj^ = l.l, 
K2=1.0, K^=-0.9, K^=0.6, A^=0.25=A2. 


In the Figs. (4.2 

and 4 

.3) we 

see the 

effect of 

simpl e 

microfluid parameters 

with 

Q=0 

(for zero 

flow) and 

0=0.5 

respectively. These 

f igures 

show 

that pressure rise 

(-Ap) 


increases as the amplitude ratio (£) increases and the increase 
is sharp for £ > 0.3. Further, the increase in parameter , and 
decrease in reduce the pressure rise, while, the parameter E^ 
has only a marginal effect on Ap. The reduction in pressure rise 
due to the variation of the parameters E^ and E^ is more in the 
case of zero flow. Pressure rise gets further reduced as the time 
averaged flux increases. These effects are elaborated in Figs. 
(4.4 to 4.6). As expected, the maximum pressure rise occurs at 
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FIG. 4.2 VARIATION OF (-AP) 
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6 FOR DIFFERENT E^ 






FIG.4.3 VARIATION OF (-AP) WITH 6 FOR DIFFERENT Ei£ 28 .E 3 . 
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FIG. 






FIG.4,6 VARIATION OF(-AP) WITH & FOR DIFFERENT Ei a E 3 . 
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zero flow rate (Fig. 4.6). Fig. (4.6) also shows that zero 
pressure difference across a wavelength is achieved around Q=0-5 
and the value of Q for which Ap=0, becones higher as 
increases. The effect of and gets reversed for Q > 0.6. 

Friction force at the wall ^ vs, £ (the anplitude ratio) is 
plotted in Figs. (4.7 and 4.8) while in Figs. (4.9 and 4.10) we 
have shown its variation with and Q respectively at ^ = 0.4. The 
effect of various parameters is similar to the effect obtained in 
the case of pressure drop, however, it is noted that the effect 
of and E^ gets reversed for Q > 0.45, e.g. , Fig. (4.9) shows 
that the friction force increases with the increase in E^ at Q=0 
but decreases at 0^=0. 5. The friction force ^ varies linearly 
with time averaged flux and has almost linear variation with 
respect to E^ for Q=0.5. 
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FIG.4.7 VARIATION OF (-5^) V^ITH 6 FOR DIFFERENT E^ & 1^ 
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FIG.4.8 VARIATION OF (-^) WITH € FOR DIFFERENT Ei & E 3 . 
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4.3 Part II s FLOW THROUGH A CHAMHEL 


In part I we have studied a siiaple iBodel for the peristaltic 
transport of a simple microfluid by considering the flow to be 
inertia free. In this part ve wish to analyse the problem for low 
Reynolds number flow. It is well known that the problem of 
peristaltic motion in general involves four parameters i.e. (i) 
amplitude ratio, (ii) slope parameter, (iii) Reynolds number and 
(iv) dimensionless time averaged flux or pressure rise per 
wavelength [Jaffrin and Shapiro (1971)]. It is not possible to 
get a general analytic solution for arbitrary values of all these 
parameters, hence, approximate analytical solutions have been 
attempted by many authors assuming one or more of these 
parameters to be small, [Burns and Parkes (1967), Barton and 
Raynor (1968), Fung and Yih (1968), Zien and Ostrach (1970), Devi 
and Devanathan (1975), Kalmal (1978), Radhakrishnamacharya 
(1982), etc]. In view of this we consider here peristaltic motion 
of a simple microfluid through a channel under the assumption 
that the wavelength (X.) of the peristaltic wave is large in 
comparison with the mean half width (h) of the channel i.e. the 
slope parameter (j&* = h/X) < 1. Thus approximate solutions for a 
stream function is attempted, using regular perturbation method 
with £* as the perturbation parameter. The analysis is valid for 


low Reynolds number. 
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4. 3> 1 Formulate ion s 

Let ua consider the moving frame of reference having 
velocity of the wave speed C and moving in the positive axial 
direction. The coordinates in the stationary wave are 

transformed to those in the moving frame of reference (X,Z) as 
follows [Fig. 4.1(a)] : 

• . 

X = X - Ct and Z = Z . 

f ic * 

Thus, the equation of the channel wall [(Z = n (X ,t)] in the 

moving frame of reference can be written as, 

Z = n (X) = h+bsin(-^^) (4.27) 

where 2 h is the constant mean width of the channel, b is the 
amplitude and X is the wavelength of the peristaltic wave. 

Now we Introduce the following non-dimensional variables. 
(u,w) = (U,U)/C, !,>..= ^ ., p = p’h/pc\ ,x = ^, z = Z/h, 

1 J L 1 J ^ 

r) = Tf 7 /h, c? = b/h, £ = h/X, R = phC/p is the flow Reynolds 

number, [where (U, 0 ,W) are the velocity components in the moving 

» 

frame of reference and v. . - the gyration tensor in the moving 

* 

coordinate system has only two non-zero components, namely ^^3 

and vjjl. 

Using the above non-dimensional scheme the equations of 
motion for flow of incompressible simple microfluid in the moving 
frame of reference can be written as follows : 
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(4.29) 
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Sx ^x Sz Sz 


Six * Sw « 

2(Ei + 2E2) + E3 2Vj^3^ " ^ ^ ^ = ° 


(4.30) 


*2 *2 ^ ^31 
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^ ^ ^ ^ Sxc Sz Hz 


- ZCEj+ZEj) * E 3 2 «j33j - « 


tfv ^ ^ _ 

^ dz 


0 (4.31) 


alonewith the equations of conservation of mass i.e. Uj^j^ 


= 0 (where E^^ , E^ , E3, 


ic. are as defined in Chapter I & II) 



Ill 


The non dimensionAl boundary conditions in stoving frame of 
reference can be written as : 

At the channel wall, we have no slip condition for velocity 
component, where as the gyration tensor is assumed to satisfy the 
condition as given by eqn. (1.25) in Chapter I. 


(i) 

u = -1 





1 



(ii) 

^(13) ^ ‘^1*^13 

Ii 

* 

Sz 

xlc 

1 




(Hi) 

‘^[13] "" ■^2'^13 

1 

_ "^2 

2 

Su * 

J 

Sv 


► Z=r)(x) 

(4.32) 




1 

- 

- 




(iv) 

n 

0 


at z = 0 




(4.33) 

(V) 

"l3 == “31 ° “ 









• * 

Further we prescribe flux (20 ) in the moving coordinates. This 
In the non-dimensional form gives, 

YjCx) 

J u dz = q (which is constant) 
o 

[Note that the flux 2Q* in the moving coordinates is related to 

* 

2Q in stationary system by Q = 0 + bCh and is constant]. 

Now introducing the stream function v* 


u 


iz 


* ^ 

and V = ^ 
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the equations (4.28 to 4.33) can be rewritten in terme of stream 
function after el iminatine P, as 


*£* I ^ ^ ^ ^ f 9 


Sz Sk I . 2 

Idz 




2 *2 "ll 

& y/ 

7? \\ 


2^ 2 
Sz Sx 


( 2 .E 3 ) 


s\ ^ 




. . f IZmi IS132 

1 Sz Sx^ 


1 r 

rn- 


“I3H iSml - n 


(4.34) 


0-2 - -2 *2 .2 
Sx Sx Sz Sz 


2 CE 1 - 2 E 3 ) .^^ 3 ^ . E3[2v^^33 =0 (4.35) 


_ _ *2 sK _ _ *2 sK _ _ sK _ _ s\ 

(K +K3)-. |i + T- * (K -K ) + CK -K ) |i 


"(31) * ^3 ^*'[31] " ^ ■“ 


The corresponding boundary conditions are 


±JL = 0 


at z = 0 


(4.37) 


"13 “ "31 ' ° 
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Sz 


= -1 


at 


2 = T) 


^"(13) = ^ 




[13] 




fV 

ax2 


at z = yj 


at z = 7) 


(4.38) 


The condition prescribing flux in the moving frame yields. 


V' = q at z = n . (4.39) 

In the following, we attempt an approximate solution of the 

above system using the perturbation method, with perturbation 

* 

parameter being the slope parameter . Thus 


y; = V'q + ^ + 0 (£ ) 

* *2 

v_ = M +£ M, + 0(c ) ► 

13 o 1 

* *2 

V, = N + N, + 0 (-? ) J 

3 1 O t 


(4.40) 


Using the above expressions in eqns. (4.34 to 4,36) and assuming 

that R*^* 0 (£*), we collect the coefficients of like powers of 
£ , which give the coupled equations governing yf , M and N of 
different orders as follows. 
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ZEROTH ORMIR 




9 (M +N ) 


( 2 +E 3 ) ^ 

Sz 


o . ^ _'-'o - o' _ , <» 


»z 


'1 ~T 


= 0 


(4.41) 






(K.-K,) — ^ +CK -K ) ^ -(E +2E )Ctt +N )+E_ 

1 ^ ^ Sz 1 2 o o 3 




(N -K )+ 


« ** 1 • 


=0 


(4.42) 






oZ Z 


(M -N )- 






= 0 


With corresponding boundary condtiiona 


(4.43) 


V^, 


/hi/ 

c 

'WT 


V', 


= q 


(13) 


[13] 


at z = 0 


at z = n 


at 2 = T) 


.2 

•i (M +N ) = A — ^ at z = T) 

2 00 1 ^^2 


.2 


(M -N ) = A, 


2 - o * "2 ^^2 


at z = Y) 


(4.44) 


FIRST ORDER 




C2*E33^ -Ei^CMi*Nj)-E3^(Mi-IIj) = 2r’ 
dz Sz Oz 



^ W 


Sz 

ifxSz^ 

dx * 3 
Sz 


(4.45) 
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(Ki + K3)-^ + -(Ej^ + 2E2)(n^+N^) + E. 






&z 


= 0 


(4.46) 


9 9 

^ M d ^ N 


a\ ■■ 

(„1-Ni) . ^ 


Sz 


(4.47) 


The corresponding boundary condtiiona are 


V. =0 




Sz 




V 


v> 


(13) 


[13] 


(M^+N^) 


(«,-N^) 


at z = 0 


at z = n 


at z = ri 


*2 

s v. 


= A 


1 dz2 

s\. 


= A 


2 ^z2 


at z = 7 ) 


at z = 7 ) 


(4.48) 


Integrating eqn. (4.41) twice and using the condition at z - 0; 

we get. 


.2 


Sz' 


(2 + E 




(4.49) 


(where A is the integration constant). 
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The eqns. (4.42) and (4.43) on using above equation can be 
written in the following form .- 


K. 


dz2 

2 

■ "1 

M 

0 

^ "2 

■ d^ 

2 

" "2 

dz^ 

2 

■ ^3 

M 

0 

- ^4 

■ d^ 

2 

■ "^4 


N = 0 
o 


N = 


AE^z 

( 2 +E 3 ) 


(4.50) 


(4.51) 


(where are as defined in Chapter III). Now equations (4.50) 
and (4.51) can be solved in the usual way as in Chapter III. Thus 


solving for M , N and y we get, 

0 0 o » ' 

^ r 1 

*0 ‘ 1 ■ ^ F.(fl) G - 

i = l *■ i 


z - 


M 


^ f ^^1 

2 

i = l I” 


Cn)z - P,(a^-a^) sinhoi z G (r?) 

^ 1 D IX 


, 3 , 2 . 

^ (z -3 y) z) 

2r, 

(4.52) 

3Pj (q+r})z 


n 


t ff 3P F (Y))z ) 1 3P 

2 •||^sinha^z - — ^ G^(r))j + — j- (q+'r ))2 


(4-53) 


(4.54) 


where Fj|_(z) 


sinha . z 

— - z cosha.D 


e . 
1 


Gj^(rj) = 3(q+Y)) j^P^(l- 2 Aj^- 2 A 2 ) + P^(l + 2Aj^-2A2)j / |^2y) (^ 3^5 ^4^5^j 

G2(r)) = -3(q+7))[P5(l-2A^-2A2)+P3(l + 2A2-2A3)] / [2Tj^(P3P^-P4P5)j 
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P 


2 + i 


(v?) 


P 


4 + i 


(17) 


^ 3*^1 ' [«V(2 + E3) (K^K^ - K2K3)] 


(K 1 K 4 - K 2 K 3 ) 



(4 

= sfi - A -A 1 

F.(n) 

- j^sinha^yjjj 


^[2 1 2 J 

2 

n 


= -aj^l + + A^jp.CT)) + [d - 


( 1 = 1 , 2 ) defined in Chapter III (Eqn. 3 . 16 ) and 


a 


2 

5 


K.K- et“ 
14 4 


(KiK, 


K 1 K 3 « 


- K 2 K 3 ) 


2 

2 


[Note : Throughout this chapter the summation index ,i, takes the 

- - 2-2 

value 1 and 2 ]. Also = ca, and ot , ft are the roots of 

the eqn. ( 2 . 38 ) . 


Using eqn. ( 4 . 52 ), we find Vq^zz’ ^ozzz’ gives : 


-2 


oz xzzz 


ainh(a^z) 


W W 
^OK^OZZZ 




) = z[q2p^+qP3+pJ + 2^pPlO^‘»^ll'^^2] 

‘*^14 ^ ^i5 ] ^ cosh(a.z)[q^P.^ + qP.^ + P.g] 


+ z^sinh(a^z) + qQj ^2 * 0 ^ 3 ] ~z^cosh(a^z) + qQ^g + Q^J 

- (- 1 )^"^^ ainh(a^z) coah(a^^j^z) [<l^ 0 £ 7 ‘*^' 5 Qi 8 ^®i 9 ] j ( 4 . 56 ) 
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where to ®ij known functions of r> and rj . 

Integrating equation (4.45) twice with respect to z and using the 
boundary conditions at z=0, we gat 


(2.E3)^ 

Sz 


(2^E^) 


CM^+N^)- 


( 2 .E 3 ) 


= B(z) 


(4.57) 


z z 


2 

R f r 

0 

S ^ 

0 

0 

^ v; 

0 

2 -^E 3 

j j 

0 0 

Sz 

SxSz^ 


az^ 


(dz)(dz )+Jz 


(4.58) 


where J is integration constant which is to be found out later. 


Equations (4.46) and (4.47) along with equations (4.57) give : 


-E^B(z) 

^ . 2 27 

K3C<»3-«i) 


K 3 K 3 



Also we get, 


(D^-oi^)(D^-/9^)N^ 


r 2 2 i 

• E3B(z) 

P -«i 


L U 

(KjKj-iCjKj) J 


(4.59) 


(4.60) 


Solving for and using the boundary conditions K^=N^ = 0 at 

z=0, we get 

M +N = i - + PaqJz + ■*■ ^*^32 ^33^ ^ 

K 3 ( 0 . 3 - 0 .^) 

+ (^^^34 + ^^35 * ^^ 36 ^ ^*^^^37 ^^^38 * ^ 39 ^ ^ 
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- Aj I * 2 


1 *- 


z cosh(a.z)Cq 


+ z ainhCa^z) (q Ui7 +qUig+U^ 5> ) + z^ cosh(a^z) Cq^X^j+qX^2^^i3^ 


+ z^ sinh(a^z) (q^X^^+qX^5+X^^) + sinh(a^+a^^pz (q^X^^+qX^g+X^^ 


+ sinh(a^-a^^^)z(q Y^^ + qY^2'*‘Yi3) 


'[- 


+ Wj- Ag(a^ - otg) + Ij sinh(az) 


[- 


+ P |- A, (/ 9 ^ - oig) + 1 j sinh (ftz) 


(4.61) 


-EgB(z) 

”r*^l = ; + P^qJz + (q P^i + qP 42 + ^43^ 

KsCot^-^i) 


+ (q P44 + qP45 + qP4^) z + (q P;,, + qP;,^ + p>,o) z 


47 ^ 48 49 


- A3 J sinhCa.z)(q^U^^+qU.2+U.3) + J 


i »- 


z coah(a^z)(q '5^i4'^9Xj. g+Y^^ ) 


+ z^ sinh(a^z) (q^^^^+qY^g+Y^ ^ ) + z coaZ(a^z) (q‘ ^il ■"‘*^ 12 ^^ 13 ^ 


+ z^ ainh(a^z) (9^2 j, ^+qZ^g+Z^^) + 8inh(a^+a^^^)z (q Z^^+qZ^g+Z^^) 


+ sinh(a^-a.^j^)z(q' AjjtqAj2 +Ai3) 




+ HI- AgCoi^ - ag) - ij sinh(oiz) 


[-* 


+ P I- A3 (/^^ - otg) - 1 j sinh iftz) 


( 4 . 62 ) 
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(K3K2 - 

where A- = , 

K 1 K 3 

M & P are constants of integration and are yet to be determined, 

while j ’ ^ij’ ^ij known functions of T) and 

its derivative. 


Substituting and from eqns. (4.61), (4.62) in eqn. (4.57) 
and integrating the resultant equation with the boundary 
conditions 


= 0 at z = r) and = 0 at z = 0 , 

we get, 

3 2 3 2 

'"1 = "^50 « J- ^ ^51 * ■’‘^52 * <-6 2- ’ 

+ Cq^PjjtqP^S + Pjj) c| ♦ ('l^P 57 *‘)P 56 *P 5 j) (|j “ 


* (a^p^, - 0P58 * P55) ( ll - V > 


2 


1 


sinh(a^ 2 ) zcosh(a^z) 


2 
a . 
1 


q + q ^12 * ^13 


2 


i »■ 


zcosh(a. z) 


a . 


2 sinh(a^z) 

3 

a. 

1 


Z7jsinh(a.r)) zcosh(a.Y 7 ) 

JL - <*• 


a. 

i 


q E.^ + q E .3 + E.^ 
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z sinh(a.z) 
1 


a , 
1 


4zcosh(a^z) 6 sinh(a.z) T)^zcosh(a^z ) 


3 

a. 

X 


4 

a. 

1 


2y^zsinh(a Y) 2zcosh(a .r? ) 

+ z ^ 

2 


3 
a , 
1 


q E +qE +E 
^ 17 ^18 19 


2 


r 3 2 

z coshCa^z) 6 z sinh(a^z) 18zcosh(a.z) 24sinhCa^z) 


2 

a. 

1 


3 
a . 

X 


4 
a . 

i 


5 
a . 

X 


3.2 

T) zslnh( a^Tr 7 ) zcosh(a^Y)) 6 'j 7 ZslnhCa .y? ) 6 zcosh(a^z) 


2 

a. 

1 


* ‘>^12 * ‘’ i3 


2 


4 3 2 

z alnh(a^z) 8 z cosh(a^z) 36z slnh(a^z) 96zcosh(a^z) 


3 

a, 

1 


4 

a, 

1 


a , 


1 2 0 s 1 nh ( a z ) 


4 3 

zr? cosh(a.y)) 4y) zainh(a yj) 

1 + S 


12r} zcosh(a^-rj) 


24rjzsinh(a.T)) 24zcoshCa^'irj) 


5 

a. 

1 


1 ®i7 * '’®ie ®i9 


] 


2 


r sinh(a +a .) cosh(a. +a^^j^ ) y> 


i •- 


1 i + 1 

(ai+a.^^)‘ 


c^i + 


^ ^14 ■" * ^i6 
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r elnh(a^-a^^ 

r (a.-a. , ) 

i *- i + 1-^ 


i + cosh(a.-a^^j)y7 


^ ^17 ^ " ^19 


„ „ 1 sinh(otz) _ z coah(oiy)) 

" 60 -2 ■ ^ 

OL a 


+ P P 


61 


sinh(/9z) _ z coshC/Ty?) 

ft 


(4.63) 


The remaining boundary conditions at z = r)Cx) determine P, H and 
J as follows : 

P = J P * q^'Pjj + qP,4 ^ P„ * P,8 I (qVjj t qA*^ + a'^) 


M = J 


80 


81 


^79 I 

2 * 

(q A.^ 

* '>'‘15 

i 



+ P.S.^ + 

83 

^84 2 

2 * 

(q A.^ 


i 


* ''85 1 

2 * 

(q A^4 

+ qA*5 


(4.64) 


(4.65) 


^ = ^86 ^AA ^ ^89> 


87 


88 89- 


+ q A*8 + A*9 


+ ^ sinh(a^r)) q^ A^. 

1 

+ ^ cosh(a^17) q^ + q ^^2 * ®i 3 ] 

[ O ic * * 1 

<3 ®i4 ®i5 * ®i6 J 
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. I coah(»^ * [ ,2 b*, * ^ he * ] 

* 1 C^j ♦ q . C*3 ] 

* 2 <4 + ■> <5 * <6 ] 

* 2 ko ^84 - !2£2S^ 

L V. ^ ® 


61 78 


f — 

sinh(/9r7) 

' -2 

•r)cosh(/ir)) ] 



ft J 



* 2 1^40 ^35 

* 1^ ot a 


sii 
1 

u 


. ^lnh(^) _ YfcoBhCfSr}) 

^61 ^90 1 ^2 ^ 


}] 


(4.66) 


STRESS AT THE WALL 


The atreaa at the wall la given by 


13 


f 2 - ^ ft33 - h.) IT 


1 ^ 


(4-67) 


From the eqn. (4.67) after omitting the second and higher order 
powers of £ , we get the stress at the wall as 


" '*33 - *11^ S- “ ^ 


(4.68) 
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Thus, we get 




= f 1 _ 

[Uz^ 


p *p r i 2 a^wl 

^1 (13) ^3l^[13] 2^2“ ^ 21 


Sz 


ax 


(4.69) 


'Sc> (’ 


13 11 


= - 4 ^ 


* *2 
* y/ 


szax 


Substitutlne these value in eqn. (4. 68) we get 


- - Zli 

w (jC 


(1 


^ 1^1 * ^ 3*’2 




+ £ 


az 


a 

1 T~2 

az 


at z = ■) 7 (x) 


(4.70) 

This value is numerically calculated for different set of para- 
meters of simple microfluid using the expressions for y' and y/. . 

o 1 

PRESSURE DROP 

Eqn. (4.20) in terms of stream function can be written as 


- ££ + ^ ^ 
ax 3 

az 


-3 

+ _1 ** 

2 




2 2 

. [ 


R -r 


az 


az azax ax 


az 




(4.71) 


which gives the equations governing zeroth and first order 
pressure as follows : 
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Sx 


11 ^ 

2 az 


^3 ^ ^1 J* 

V o IS 


[^-"o] * ” * 2^> 73^ - -r i; 


^z" 


(4.72) 


3^ 


E.. S'y^ E. 


5i ["i-nJ + (1 * ^ ^ I- (BjtNj) 


2 ' * ^ 
Sz 


[ V' V 

oz ^c 


- R iv V - V V 

oxz ox ozz 


Therefore, we get 


(4.73) 


££ = l!£ » f!i 

Sx ~ Sx ^ Sx 


Cr?) 




P - — 
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(G.) eTG- 
1 11 


a , 
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^2 2., 

- G . © . 1 

+ 1 _ '9 _i_£ 

7 A 2 
r, . ot. 


+ q 


* 

£ P 
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3 fv 1 * *r • 

% \l -i-3-i - 1 ^ ^ R ^ (T)) (P.-Zr^F^P^.SP^P^r,) 

>) li 1) J I 


r , G. I 

-2 -4 

l rj J 


2e‘G.(r)) 


5 2 
n a. 


I 


3e^G^(Y)) 
1 1 


6 2 
1 i >) *i 


- (D) 


•■’■I ] 


* 

£■ P 


93 3 

r) 
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t * 

+ I? R 


f 


(17) P, 


PiP^tXtd) 


+ 2P2P3i7Cr)) 


1 - 


(Gj) 


'■ l“‘] "77" ■ [ “■’ 

■'•ll 




2^2 

ri ft 


27) (77) 


+ 0 (£ ) (4-74) 
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The pressure drop across the wavelength thus, can be evaluated by 

1 

" " j ^ If ^ (4.75) 

o 

Using eqn«(4,75) Ap is also ntimerically calculated and its values 
are plotted for different values of the simple aicrofluid para- 
meters , 

4. RESULTS AMD DISCUSSION 

Ue have presented approximate solution for a stream function 
as an asymptotic senes in terms of slope parameter . It may be 
pointed out that the zeroth order solution gives the limiting 

t 

solution in the limit ^ -► 0 which corresponds to the situation 

of very long wavelength. Further in order to discuss the results 
qualitatively the expressions for pressure rise(-Ap) and stress 
at the wall i.e. have been calculated numerically taking 

K^=-0.9, K^=0.6, A^=0.25, A2=0.25 and q=0.2. The 
results are plotted in Figs. (4.11 to 4.20) for various values of 
di f f er ent parameters . 

In Figs. (4.11 to 4.15) the effect of the simple microfluid 

parameter and the Reynolds number is shown. As in the case of 

part I, pressure rise increases with the amplitude ratio and this 

* 

increase becomes very sharp after ^=0.4,, when R -1.0. Also 

preaaure riae increaaea aa the parameter increaaea, and it 

reducea with the increase in the parameters and E^ . 


The 
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FIG.4.12 VARIATION OF (-AP) WITH € FOR DIFFERENT E,,E2&E3. 
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4.13 VARIATION OF (-AP) WITH FOR DIFFERENT E2,E3 a R"^ 
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FIG.4.14 VARIATION OF (-AP) WITH E2 F’OR DIFFERENT Ej.E3&R. 
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FIG. 4.15 VARIATION OF(-AP) WITH E 3 FOR DIFFERENT Ei,E2&R. 
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variation with the aimple microfluid parameter becomes very small 
at R =1.0. These effects are further elaborated in Figs. (4.3 to 
4.5) which are drawn for 4>=0.A. The variation of Ap with is 

more, compared to its variation with respect to E. and E. . The 
effect of the parameter at R*=1.0, depends upon the values of 
Ej and E 2 o-fi- pressure rise gets reduced with the increase in E^ 
for E 2 < 0.2. 

In Figs. (4.16 to 4.20) we have drawn shear stress at the 

wall across a wavelength (for c=0.1). The shear stress profile 

changes singi f icantly with the variation in E, . Ue also note that 

the shear stress at the wall increases as the amplitude ratio and 

the slope parameter increase. The variation is more significant 

for R =5; it is observed that r increases with increase in 

w 

parameter at R =0 & 5, however the effect of the simple 

microfluid parameter E 2 on is affected by the Reynolds number. 
* 

At R =0, T shows a marginal decrease with an increase in E . But 
w ^ 

it shows significant increase with an increase in E^ at R =5. 
Thus the Reynolds number plays an important role here. 


4.. 4. CONCLUSION 

It is noted that the pressure rise, friction force and the 
stress at the wall increase with increase in the amplitude ratio 
of the wave. Also they decrease in magnitude while the simple 
fflicrofluid parameters E^ and E 2 increase but increase with 
increase in the rotation parameter E^. It is also noted that the 
pressure rise and the friction force decrease as the flow rate is 
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X - COORDINATE > 

FIG.4.16 VARIATION OF t^ALONG THE AXIS FOR DIFFERENT 
El ,E2 S' ^3 • 
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FIG.4.18 VARIATION OF rwVHITH € FOR DIFFERENT E, .E 2 8 . E 3 . 






FIG.4.19 VARIATION OF TwWITH 6* FOR DIFFERENT E,,E2 & E3 . 
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FIG.4.20 variation OF '^WITH €*F0R DIFFERENT E,,E 2 & E3 
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t - 

increased. At R -0, shows only a marginal decrease with an 
increase in but it significantly increases with increase in 

Tk 

at R =5 . Thus we can conclude that the Reynolds nuiiber plays an 
Ifiaportant role in the peristaltic motion of a simple aicrofluid- 



CHAPTER V 


SELF PROPULSION OF SPERMATOZOA IN MICROCONTINUA . EFFECT OF 
TRANSVERSE WAVE MOTION OF CHANNEL WALLS 

B. 1 I MTRODUCTI C«4 

The Spermatozoa are tiny living cells having their own 
motility. They move towards the oviduct through the mucus 
filling the cervical canal by sending waves of lateral 
displacement down a thin tail or flagellum. In recent years many 
research workers have analytically investigated swimming of micro 
organism in liquids. Taylor (1951) in his pioneering work on 
locomotion of an organism, possessing a single flagellum (e.g. 
spermatozoa) , modelled the micro organim as a two dimensional 
infinite extensible sheet of zero thickness. He thus considered 
the motion of fluid near the sheet down which the waves of 
lateral displacements are propagated. Hancock (1953) studied the 
propulsion of a thin circular filament in infinite viscous fluid 
to obtain the velocity of propulsion of spermatozoa. Reynolds 
(1965) considered a two dimensional model and used Taylor's sheet 
approach to consider the effect of rigid walls adjacent to it. 
Brokaw (1972), Shack & Lardner (1974), Shack et al . (1974), Katz 
(1974), Lighthill (1976), Shack et al . (1978, 1988) and others 

have also attempted to study the motion of spermatozoa in 


cervical canal. 
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In all these works the problem has been treated from 
continuum point of view. However, there exists a strong 
likelihood that the cervical mucus be considered as a micro 
continuum fluid. Odebald (1959, 1962) has observed that the 
cervical mucus la a suspension of high molecular weight macro 
molecules in a mesh, having a spacing of 0.03 ^jm in them, Elstein 
(1971), Davajan et al . (1971). This motivated Sinha et al . 
(1982) to consider the micro continuum approach for studying 
self propulsion of spermatozoa in a micropolar fluid filled in a 
rigid channel. In a recent study, Shukla et al . (1988) have 
modelled the cervical canal as a channel with distensible walls. 
They have considered transverse waves of finite amplitude 
travelling down the flexible walls to account for the mascular 
activity. 

In view of this we discuss here the self propulsion of 
spermatozoa through a microcontinuum filling the channel, with 
flexible walls. This chapter consists of two parts, in part I 
mucus has been modelled as micropolar fluid model while in part 
II, simple microfluid model in considered. 

5. 2 FORMULATION 

Consider the motion of spermatozoa in cervical canal whose 
physical configuration is shown in fig. (5.1). The cervical 
canal is modelled as a two dimensional channel with flexible 
walls. It la assumed that the sinusoidal waves of finite 
amplitude are imposed along the flexible walls of the channel in 
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FIG.51 PROPAGATION OF AN ELASTIC SHEET 

SWIMMING THROUGH A TWO DIMENSIONAL 
CHANNEL WITH PERISTALTIC MOTION OF 
WALLS 
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the direction 

opposite to 

the 

motion 

of 

spermatozoa. 

The 

spermatozoa is 

assumed to be 

an 

infinite 

flexible 

sheet 

which 

propels itself 

forward by passing 

sinusoidal 

lateral 

waves 

down 


its length in the backward direction. The sheet is considered to 

be swimming with the propulsive velocity V in the negative axial 

P 

direction. Further we assume that the waves travelling along the 
channel walls and that along the sheet are in synchronisat ion 
under the steady state (Shukla et al . 1988). Thus they have same 
wave speed C (along positive axial direction) and wavelength 

f f f 

In a fixed frame of reference (X ,Z ,t ) the forms of the channel 

r 

boundaries and the sheet at an instance t are given by 


r t t 9 r 9 9 f 

Z = h + b, sin ^ (X - Ct + V t ) = H, (X , t ) (S.la) 

IX pi 

Z = -h + b^ sin (X - Ct + t ) = H.(X ,t ) (5. lb) 

z’ = bp sin ^ Cx' - Ct' + Vp t’) = Hq Cx',t’) C5.1C) 


where the amplitudes of the sinusoidal waves travelling 
along the upper wall and the lower wall respectively, bg is the 
amplitude of the lateral wave along the sheet and 2h is the width 
of the channel. 


It may be noted 

that the relative velocity of 

the waves 

with 

respect to the fixed 

» 

frame is Hence in a 

moving 

frame 

f » f * 

(x ,z ,t ) where [x 

=x'-Ct ~ Z], 

moving 

with 

9 

the velocity 

the positive axial direction. 

the sheet 

and 
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the walle appear to be stationary and the flow is steady. Thus in 
this moving frame of reference, the channel walla are given by 2 

f 

= h.h^(x), z = h.h^Cx) and the sheet by z = h.h^Cx), where 

hiCx) = 1 + sin 2nx ( 5 . 2a) 

h 2 (x) = - 1 + ^ sin 2frx 

h^Cx) = £ sin 2rTX 
u o 

b. 

^ i ~ H — * ^ ~0,1,2, x = 

Further the Reynolds number of 
the sperm, density and viscosity of 
velocity of propulsion is found to 
studies. Therefore we consider 
predominately due to viscous forces 
been neglected. In the following the superscript ( + ) refers to 

ft y 

various quantities in the region Hp s: Z i while (-) refers 

f ft 

to those in the region £ Z ^ 


C5.2b3 

(5. 2c) 

X /X . 

motion based on the length of 
the surrounding fluid and 
be of order 10 ^ in such 
here that the motion is 
and the inertial forces have 


B. 3 PART I I FLOW THROUGH MICR<M'Ctf-AR FLUID 
S. 3. 1 Governing Equations 

Here we consider mucus as a micropolar fluid. The equations 
governing the inertial free flow under the very long wavelength 
approximation have been deduced by Sinha et al . (1983) and are 
smimilar to those in lubrication theory for micropolar fluid 
[Prakash and Sinha 1975]. Ue reproduce these equations here in 
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the non- d i inens i onal form with respect to the moving frame of 
reference. Thus for the region h^Cx) z ^ h^(x) we write the 
equations as 


. (1 * ^ * s III 


( 5 . 3 ) 


££- 

Sz 


= 0 


( 5 . 4 ) 


2 ,+ 


■ " ■' "" ■ y 


- E 


Sz 


2lc ^ 


= 0 


( 5 . 5 ) 


where the following non-dimensional scheme has been used, 


X = X /X ; z = z /h ; t = Ct /X ; V = V /C 


4- ^ / rn 

U = U /C 


jl K/r. _ P 


( 5 . 6 ) 


V = 


V 


tJh‘ 


xfy^ 


’ + ’ + ’ + 

[Here the p is the pressure, u the axial velocity, ^ the 

inicrorotat ion vector ^ are viscosity coefficients for 

inicropolar fluid]. 

It may be noted that ic is coupling parameter as the equations 
become uncoupled as k-»0 and v is material coefficient parameter. 
The equations governing the flow in the region h 2 (x) ^ z h^fx) 

can also be written similarly. 

The boundary conditions with respect to moving frame of 


reference are 
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u = u 

= -1 , 

+ 

4 ’ ~ 4 ‘ 

= 0 at z = 

^0 

(5.7) 

= V 

P 

- 1 , 

= 0 

at z = 

^1 

(5.8) 

u“ = V 

P 

- 1 . 

= 0 

at z = ^2 


(5.9) 

The boundary conditions 

for , 

^ denote that 

the fluid - 

solid 


interface interaction is so strong that the microstructur e (fluid 
particle) does not rotate relative to the wall. 

Further the sheet is self-propelling and the forces exerted 
by the fluid on it must balance for its motion to be with a 
constant velocity i.e. 

J (T^ + T") dS = 0 (5.10) 

S 

where and T~ are the resultant of forces acting on the upper 
and lower surfaces of the sheet respectively and S is the surface 
area of the micro-organism [Shack & Lardener (1974), Shukla et 
al. (1988)]. 

S. 3. AnalyslB 

Solving the equations (5.3) and (5.5) using the boundary 
conditions for we get 

1 

J z 




= P 


hj^ cosh a(z-hQ) ~ hg cosh <»(z-hj^) 


sinh otH. 
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^ ^ Ch^ . - (-1 


‘"“J 

f 

cosh oi(z“hp) " cosh ot(Z"-h^) 


2 

— 



1 



L J 



(5.11) 


+ V 


cohh aCz-hjj) - cosh 01(2 ” 

Fj^(hj^) sinh oiH^ ^ 


Substitutinfl for 4> from eqn. (5.11) in eqn. (5.4) w© get 


+ r> + 

u = P 


7 (z^ - zhp - 2h^ + h^hp) - 


Hi « 


2L sinh otH, 


cosh oiCz-hg) 


+ cosh oi(z-h^) ~ cosh otH^ - 1 


F.Cz) 

+ V — - 1 

^ F ^( h ^) 


(5.12) 


222 k 4 

where = N L , N = (_51_) , L = (-^) 

2+k u 


1/2 


(5.13) 


FjCz) = C*-h„)- ^ 


L sinh oiH- 


-Icosh otCz-h^ )-cosh oi(z-h^ ) + cosh 


( 5 . 14 ) 


Hi (hj^ - hg) and P - 


Procesdinii in ths similar way for ths othsr rsfiion also, i.s. 

(h«(x) z h (x)) we obtain, 
z o 
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u = P 


7 Cz^ - zhp - zhg + h^hp) - ^ 




2L sinh aH, 


cosh otCz-hg) 


+ cosh c*(z-h^) - cosh - 1 


F,(z) 

V — I - 1 (5.15) 


= C!:-hu)- -j— 


L sinh otH, 


cosh oi(z-hQ)-cosh ot(z“h 2 )‘^cosh 


(5.16) 


“2 = (^z - “■* ='■ - ^ 


The dimensionless flow flux q (= q in the moving frame of 
reference is obtained from : 


■ = I"* 


q = I u dz and q 
h. 


= 1 " 


dz 


which on usinfl eqns . (5.12) and (5.15) gives 


P G, + G, - H. 
1 P ^ 1 


P G- - G. + H 
3 P ’ * 


(5.17) 

(5.18) 


where 


oiH 


H^/12 + 


2L sinh 




2ainh 



om.. 


- cosh “ 1 
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Eqns. (5.17) and (5.18) can be rewritten and the pressure 
firadients can be obtained as. 


fp 

Sx 



(5.19) 


(5.20) 


Intefiratinfl the equation of continuity it can be noted that the 
two fluxes, q"" and q' are constants. Further, it is evident that 
Ap. the pressure rise over a wavelength is same for the two 
regions. Therefore, integrating eqns. (5.19) and (5.20) over a 
wavelength we get the following relations among the three unknown 

quantities V , q and q . 

P 


V I,, - I.o = 


- ^ hi Vl2 ■ "13 


(5.21) 



149 


^ ^21 * ~ ^23 * 


where 


11 


= J 


dx 


12 


J 


XT ’ 

1 


13 


= 1 


^ dx ; 


(5.22) 


21 


= J 


dx 


22 


' J 




dx ; 


23 


= J 


“2 


dx 


(5.23) 


Using the stress-strain relationship and long wavelength 
approximation we write the force equilibrium condition, eqn. 
(5.10) as. 

1 



(5.24) 



(5.25) 


which gives the third necessary relationship to determine V^, q"'^ 
and q“. Here [x] indicates the difference in the quantity x above 
and below the sheet. 
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Eqh . C5.25) on substituting for , gives 

Sz 


Q I,, + q I-, + V I,, + I,. = 0 
31 32 p 33 34 


(5.26) 


where 


•31 


= 1 


^ (h^ + H^/2) dx 


‘32 


= J ^ '"o 


+ H 2 / 2 ) dx 


■33 


I ['*'« 

o t- 


H, G. 

2 . 4 


7-^ U; 


»2 


+ — ) - Ch„ + “) 


o -2“^ 117 


^ ^ 
I 


7 T 7 H 7 J 


dx 


34 


■|[ 


®i h ^2 ®2 


dx 


(5.27) 


Solving eqns. (5.21), (5.22) and (5.26) we get. 


V = 


«- '•11 ^21-’ 


^ 31^13 ^ 32^23 _ 

~T 7 r ~ ■■■'i:: "" " ^34 


11 


■21 


[ 


(5.28) 


^12^3l''^ll " ^32^22^^21 ■*■ ^33 


The total flux q 


+ — . 
q q IS 


q = 




Ap + 




(5.29) 


Now, we can obtain the flux Q in the stationary frame and Q is 
related with flux, q in the moving frame by the following 
relation [Shukla et al . 1988], 
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0 = q + [C1 - Vp) (H^ - H^)/h j (5.30) 

Averaging this Q over a time period, we get the time averaged 

flow 0 terms of a and V as 

P 

0 = q + 2 (1 - Vp) (5.31) 

5.3.3. RESULTS AND DISCUSSION 

In our study here the additional parameters, as compared to 
the study of shack and Lardner (1974), are the micropolar 
parameters N and L & the amplitude parameters and ^2’ 

parameters N and L arise as a consequence of rotational viscosity 

and the material coefficient. As the parameter N -*■ 0 i.e. x: 

(rotational viscosity) -♦ 0 the results reduce to those for the 
Newtonian fluid. The limiting case of L -♦ a> also gives Newtonian 
results. Thus, the present analysis in these limiting cases, 
yields the results of Shukla et al. (1988), for the particular 
case of no longitudinal motion of the walls. The parameters 
and £ represent the amplitudes of the waves on the channel 
walls. The case of £. = is equivalent to imposing peristaltic 

motion on the walls. Uhile ° no wave motion on 

the wall. Thus, for no wave motion on the wall, and Ap=0 the 
present study reduces to that of Sinha et al . (1982). 

Ue have numerically computed the values of propulsive 

velocitv V and the time average flux Q for different values of 
P 



152 


micropolar parameters and Ap = 0.0, -0.05, = 0.0, 0.35 and 

X « 

= 0.0, ""0.25. The results are plotted in the figures (5.2 to 

5.143- The propulsive velocity ,V , versus £: is plotted in 

p o 

fifiures (5.2 to 5.5), while its variation with respect to L is 
shown in figures (5.6 to 5.9). Figures (5.2 and 5.3) are plotted 
for^^ =r 0 = 4?2 ^nd Ap -0.05. They show that for a given 

micropolar fluid the propulsive velocity increases as magnitude 
of increases. This is similar to the behaviour as observed by 

Shack & Lardner (1974). Further, at a given amplitude , the 

o 

propulsive velocity shows, significant quantitative increase 

with the parameter N. In figures (5.4 and 5.5) we have shown the 

effect of wall amplitude variation for a given micropolar fluid. 

The propulsive velocity decreases as ^ increases, and for given 

it further reduces as magnitude of increases. Thus the 

effect of transverse wave of the walls is to reduce the speed of 

the spermatozoa. In fact for certain combinations of the wave 

amplitudes the direction of the propulsion can even be reversed. 

From figures (5.6 to 5.8) it is observed that for ^ = 0 . 4 , ¥ 

o p 

initially increases with L and then shows a decreasing trend, 
i.e. the propulsive velocity curve reaches maximum and then tends 
to the Newtonian value on L oo. The value of L at which the 
maximum occurs very much depends upon the parameter N and the 
curve approaches maximum at smaller values of L as N 1. 

Further, it is observed from Figs. (5.5 and 5.6) that the effect 
of peristaltic wall motion is to reduce considerably . In the 
case of no wave motion on the walls, shows slight decrease as 
Ap decreases from 0.0 to -0.05, and this decrease becomes more 
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FIG.5.6 VARIATION OF Vp WITH L FOR DIFFERENT N. 
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FIG.5.7 VARIATION OF Vp V\^ITH L FOR DIFFERENT N. 



159 



FIG.5.8 VARIATION OF Vp WITH L FOR DIFFERENT N. 
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pronounced as L increases or N decreases (Figs. 5.4 to 5.5). Fifi. 

(5.9) shows the propulsion of a rigid sheet U =0) under the 

o 

peristaltic waves on the channel walls. It shows the opposite 

behaviour as observed in the case of £ =04 

o 

In fi£ures (5.10 to 5.14) tiine avera£e flux (Q) is plotted 

versus different parameters. The trend of Q vith respect to 

various parameters is opposite to that of the propulsive velocity 

. It is noted that Q decreases vith the increase in the 

roafinitude of and the effect of peristaltic wave on the wall is 

to increase Q considerably. The variation of Q vith the 

Hiicropolar parameters N and L depends upon the other paramet ers . 

For example, with peristaltic wave on the walls Q decreases with 

N at = 0.4 while it increases vith N when dC = 0.0. It may be 
o o 

pointed out here that the results are in comfirmity with the 
results of Shack & Lardner (1974), Sinha et al . (1983) and Shukla 
et al. (1988). 


S. 4 PART II t FLOW THROUGH SIMPLE MICROFLUID 

Sm 4« 1 Goverrriing erquations 

In this section we consider simple microfluid model for 
mucus. Thus, vith the follovin£ non~dimensionalisation scheme. 






X = 


z 

IT 


t = 


Ct 


u 


"^( 13 ) " r "^(13)^ ^ 131 } " i: ^[31] 











N 



0.4 



1 


i 
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K. 


1 




E 


1 




and = — 

3 


the equations of motion of simple microfluid (eqns. 1.9 to 1.12) 
can be written in the non-dimensional form under long wavelength 
approximation Ch << X) as follows : 


Sp 

Sx 


*2 ± 
S u 

Sz 


S V 


- E 


(13) 


1 dz 


+ E. 


, *2± 

[ 31 3 ^ 1 d u 

I 7^ 

Sz 


= 0 (5.33) 


_ 


Sz 


= 0 (5.34) 


(K1-K3) ^ (v.^^) + (K,-K,) 


2 “4^ ^*^13^ 2(E^ + 2E2) 


2v 


± 0\i 

[13] 5F 


(5.35) 


2 2 

(-31’ * ^ '“135 - ^133 


”'l*'=3> n -31 


+ E. 


_ ± SvT 

^*^[13] :5F 


= 0 (5.36) 


where superscripts (±) denote the flow quantities in the region 

h (x) < z < h, (x) and h_(x) ^ z < h (x) respectively, 
o 12 o 
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Th© boundary conditions are written in diwensionl ess font as 


+ 


u = 

u = 1 

at z = h (x) 

■*1 





o 




+ 

u = 

V , 

at z - h^(x) 





p-1 





u . = 

^-1 

at z = h^Cx) 



(5.37) 







11 

0 = V3^ 

at z = h^(x). 

h^(x) 



'^13 = 

° = ^31 

at z = h (x), 
o 

h^fx). 



The conditions of v. 

13 

are perfect adherence 

( conditions 

can b#' 

taken as the 

^ result of 

strong interaction 

of fluid and 

solid 


interface* 

Further, using the str ess-strain relationship and the long 
wavelength approximation the force equilibrium condition eq . 
(5.10) can be written as 

i[g]- ■ • 

o 



5. 4. £ Analysis 

Solving equations (5.34) to (5.36) for * ^31 “ alongwith 

appropriate boundary conditions (5.37) we get. 
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C13) 


2 ^ 


B-G. (a., a.^,) + P (a.,a.^.) cosh(a.z) 
ili 1+1 3^11+1 ^1 


B^G 2 Ca^,a.^^) + P G^Ca^.a^^^) sinhCa^z) . (5.39) 


[13] 


B^ - zP 


C2 + E3) (K30.;- 


2 .r 

1^1 

i = l '“L 


^1*^5 ^*i’®i + l^ + P G^ (a^.a^^j) cosha^z 


^ S ^®i'^i + l^ sinha^z. (5.40) 


u = -1 + Vp Cz-h^) + D 3 (z) 


•V “I r f 

■ ^ ^2 - 
J J L ^ 


(z-h^) + D 3 (z) 


2 v2 

z -h 

o 


+ D^(z); 


(5.41) 


where 


^ = 


1 o 31' 


r 

“2 = ”1 
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= ITT. 


1 


X *- 




) L-Ca.z) + G,(a.,a.^,) L„Ca.z) 




} 




= ih-. 


1 ^ 








j + 1 


1 + F3Ca.) 


FcC j >3-. »a. 3 L.Ca.h.) 

5 1 i+1 j 1 1 


Gy,^ -Ca. 
4 + 1 


a 3 = 
i + l'^ 2 


j + 1 


{ 


i - 






Fc ( j > a . , a . , 3 L.(a.h,3 
5 1 1+1 j 1 I ' 


CJ = 1 . 2 , a & 43 


Lj^Ca^z) = sinh(a^z3 - sinh(a^h^3 


L-Ca.z3 = coshCa.z3 “ cosh(a.h 3 
Z i 1 i o 


L*(a.z3 = h sinh(a.z3 - 2 sinh(a.h ) 
3 X O 1 1 o 


L,(a.z3 = h cosh(a.z3 - z cosh(a.h 3 

4 jL o X xo 
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F,(a.) 


- [l * F3<a,)] / [ 


^2 n 2 
K-oi- - K,«, 
3 3 4 4 




2 E3 F,Ca.) 






^ 5 ^^i+l’^i^ j is odd 


^ 5 ^* 1 ’ * 1 + 1 ^ j is even 


_ - - - 2-2 

“ *3 - **i ^2 ~ ^ ^ ^ positive roots of the 
eqn. (2.38). is integration constant which is evaluated 
through the application of the boundary condition 


V -1 at z=h,. 
P 1 


The flux in the upper half plane is then given by 


-- 


r h-K - 11 

+ 

“ ^ 

^ = vp 



+ P 

^ 2 - 
_ ^ 


^2 k 2 


h^Chi-h^)4D3 


I 


h^h^ 

1 o . o 


Z — ■*■ T" 


- Ch^-hp 


(5.42) 


PlVp + + P3 


(say) 


(5.43) 



172 


where 


D. 


2 

ITT. 


- ^3 


l + l5 




a.^l) SCa.) . G^(a..a.^^) C(ap 


and D 


4 


2 

Tns. 


3 ^ 

1 


G_(a.,a. ,) S(a.) + G,(a.,a.j^,) CCa.} 
3 i’ i+l' ^ 4^ i' 1+1 ^ x' 


E, ^Gy(a^,a^^^) SCa.) + GoCa.,a._,,) CCa.) 


*“i + l 


SCa^) = L2Ca.hj^) - a^Ch^^-h^) sinhCa^h^) 


CCaj^) = LjCa^hj^) - a.Ch^^-h^) sinhCa^h^) 
Similarly we can write the expression for u as 


_ 


1 



f 1 

u = -1 + V 

P 

^5- 

Cz-h^) * 


-»■ P” 

4 Cz-h^) + Dy Cz) K 


z2-h2 

-2-^ + DgCz) 


C5-44) 


where 


TiV^rrrnJTir^ 


^6 = ^5 


" ^2 ■ ^o 

2 ~ 
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> 7 ^^^ = 5 

i 


a.^l) L^(a.z) + G 2 (a.,a.^^) L^Ca.zy. 


jfCgCai.ai^l) Lj(a.z) + G^(a..a.^^) L2Ca. 




1 








G.(]a.,a- ^) — 

j 1 1 + 1 


_ C-dJ"*^ f 
zir: ■ 
i I 




® 4 +jt^i-^i+l^ = 


-T^{ 


1 - F3Ca.)l F^Cj.a-.a.^^) L.Ca.h^) 


P* ^ 

^ 5 '^^i'^i+i^ = 7 +E^ TpirrY^r^:^. 


i+1^ 2 o- 


FgC j,a^,a^^^) = 


^5^*i+l’*i^ j is odd 


^5^*i’*i+l^ j is even 


which gives flux q as follows. 


«■ = -%kl ^ - >'<,'‘'2-'‘o5*®7l] * - •'o'‘'2-»»^"“7| 
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h^-h^ 

-3- ^ ^ ^ ^8 


+ Ch„-h^) 

4 O 


(5.45) 


"4''p * ‘■s H- * 


(5.46) 


vhere 


2 Y ^ 
T^l I iT 
i 


G,(a..a.^l)S (a.) . G2(a..a.^^)C 


E, ^GgCa. ,a^^^)S (a,) + G^ (a, . a , )C (a 


6^“i '“i+1- 




2 

!+¥: 


T — E Ig (a. 
Z a . 1 3 1 

i ^ 


.ai^l)S (a.) + G^(a.,a.^^)C (a.)l 


- E. 


{6,(a.,a. 


ial>S (»i3 * '^5 


S (a ) = L Ca.h,) - a,(h,-h ) sinh(a.h ) 

1 120 lO 


C (a.) = L.Ca.h_) - a.Ch,-h ) sinhCa.h ) 

1 11* 120 lO 


From equations (5.43) and (5,46) we £et 


ap 

ax 


(a 


P3)/P2 


(5-47) 


Sx 


(q 





(5.48) 
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Also we know that 

j dX . AP = j 


Sx 


dx 


i.e. the pressure rise over a wavelength is sane for the two 
regions. We get from equations (5.47) and (5.48): 


AP = a 1^2 + 1^3 


(5.49) 


^ ^21 \ ^22 ^23 


(5.50) 


where 




dx 


:i = j 


dx 


12 


= -J 




dx 


22 


-I 


dx 


(5.51) 


13 


j pj- '** ' ‘21 

o 


-J 


dx 


Now considering force equilibrium condition (5.38) and 


substituting for and we get 

dx dx * 


<■ * "32 ■> ^ ’33 % * ’34 = » 


(5.52) 
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where 


^31 


= 1 


C ^ ) dx 
2 


P P P P 

, , 1^8 9 4 . . 

^33 " ^^7 TT" " “Fr~^ 


= 1 


32 




C Tj- ) dx 
5 


■34 


J 


P P P 

' -Tf * ^ '** 


P, = (D^-Dj) ♦ D 


'>3'"^ \ - “y'"’ 


1 h [ ’>3'-^ ] - 

]} * ^ [ “ 4 '“^ ] 


* Iz [ ]} * ^ [ ] 


^ = ■>6 ^ 


Recalling that the time averaged flux Q is given by. 


Q = q + 2 Cl - Vp) . 


(5.53) 


eliminating q''’ , q~ and q = q'*’+q” from equations (5.52) and (5.53) 
we get the expression for in terms of Ap. These quantities 
have been numerically calculated. 


5.5 RESULTS AND DISCUSS! CW 

The parameters involved in this analysis are : the simple 
microfluid parameter E, , E., E_, the wave amplitudes and 

and the pressure difference across the wavelength (Ap) . Here the 

values of V and 6 have been calculated numerically for K,=l.l, 

P * 
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K 2 = 1.0, K^ = -'0.9 a? K^ = 0.6 and ar© plotted versus in figures 

(5.15 to 5.21) and (5.22 to 5.24) respectively. Figures (5.15 to 

5.17) ©how the effect of ©iiiple liicrofluid paraiseters when there 

is no transverse waves travelling along the channel walls. As in 

the previous case it is observed that V increases with the 

P 

magnitude of a© also with Ap. In case of sero pressure 

difference across the wavelength, i.e. Ap = 0 . 0 , increases 

with the decrease in E. , E« and increase in E«. However, when Ap 

X Z 3 

= -0.05 though it shows similar behaviour for higher values for 
reverse trend is observed for small It may be further 

remarked that unlike in the other chapters effect of E^, 
variation is significant here. 

Figures (5. .18 to 5.21) show the effect of transverse motion 
of the channel walls on . Even in the presence of transverse 
wave motion, V decreases as Ap decreases, and the effect of 

p 

simple microfluid parameters remains same as has been observed 

earlier (^^ = 0 = (Figs. 5.18 a 5.19). Moreover variation of 

versus loses its symmetrical behaviour for non-symmetrical 

wave motion on the wall (Figs. 5.20 and 5.21). For given ^^=0.35 

amd Ap = -0.05, V increases considerably as increases from 

P 2 

0.0 to 0.25. 

The figures (5.22 to 5.24) show that Q increases as E^ 

increases but it decreases with an increase in and 

parameters. Further, for given = -0.25, = 0.0, 0 decreases 

Z o 

as increases from 0.0 to 0.35 and a reverse trend is observed 
1 

at ^ =0.25. In the case of = 0.35 and je 

o 1 


-0.25, Q 
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-0.25-0,2 -0.15 -0.1-0.05 0 0.05 0.1 0.15 0.2 0.25 


to ^ 

FIG.5.15 VARIATION OF VpWITH Fq FOR DIFFERENT AP & Ei . 
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FIG.5.16 VARIATION OF Vp WITH Cg FOR DIFFERENT AP & E3. 
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F16.5.17 VARIATION OF Vp WITH Eq FOR DIFFERENT AP & Ej . 



• • 


1 



- 0 . 25 - 0.2 - 0,15 - 0.1 - 0.05 0 0.05 0.1 0.15 0.2 0.25 


^0 

FIG.5.18 VARIATION OF Vp WITH Fq FOR DIFFERENT E,,E2&E3 




^0 ^ 

FIG.5.19 VARIATION OF Vp VlfITH Gq FOR DIFFERENT E,,E 2& E3. 


E2 = 0.2l 

£2=0.4] El =0.9, £3=0.3 


£1 = 0.35 
62 = 0.0 
AP =-0.05 


E2 = 0.4l 

£2=0.2] = 

5=0.2 'El = 0.3, £3=0.1 




-0.25-0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2 0.25 


FIG.5.20 VARIATION OF Vp WITH £5 FOR DIFFERENT £,, £ 2 &E 3 
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FIG. 5.21 VARIATION OF Vp WITH €q FOR DIFFERENT E^,E 2& E3. 



- 0 . 25 - 0.2 -015 - 0.1 - 0.05 0 0.05 0.1 0.15 0.2 0 . 2 ! 


^0 ^ 

FIG.5.22 VARIATION OF Q WITH Cq FOR DIFFERENT E 1 .E 2 & E 3 
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^0 ^ 

FI6.5.23 VARIATION OF Q WITH €q FOR DIFFERENT £^,£2 & E 3 . 
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FIG.5.24 VARIATION OF 5 WITH Cq FOR DIFFERENT Ei,E2 8.E3 
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initially decrease with and then shows sharp increase, while 
in the case of = 0.35 and = 0.. 25 it shows sharp decrease -as 
increases froiii-O.ZS to 0.15 and then show a sli£ht increase. 

S.6 COMCLUSIOM 

In the present study, a microcontinuum approach has been 
taken to model the mucus as micropolar fluid (in part I) and as 
simple raicrofluid (in part II). It is observed that the 
propulsive velocity increases as the departure from Newtonian 
theory increases. For simple microfluid model, increases with 
the decrease in microstr etchin£ parameters and increase in 
rotation parameter. In both the cases propulsive velocity 
increases with amplitude of the wave on the sheet. Peristaltic 
wave on the wall of channel decreases the propulsive velocity. 
This sufifiests that if by some biochemical means, a peristaltic 
wave can be imposed on the walls of the cervical canal, then the 
propulsive velocity of the spermatozoa can be reduced 
considerably, so that it will not reach the point of 
fertilization. This observation is important from the point of 


view of fertility control. 
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